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ABSTRACT 

A variety of host galaxy (bulge) parameters are examined in order to de- 
termine their predictive power in ascertaining the masses of the supermassive 
black holes (SMBH) at the centers of the galaxies. Based on a sample of 23 
nearby galaxies, comprised of both elliptical galaxies and spiral/lenticular bulges, 
we identify a strong correlation between the bulge gravitational binding energy 
(Eg), as traced by the stellar light profile, and the SMBH mass (M,), such that 
M, oc Eg'^. The scatter about the relationship indicates that this is as strong a 
predictor of M, as the velocity dispersion (a), for the elliptical galaxy subsam- 
ple. Improved mass-to-light ratios, obtained with IFU spectroscopy and I-band 
photometry by the SAURON group, were used for those sample galaxies where 
available, resulting in an energy predictor with the same slope, but with reduced 
scatter. Alternative M, predictors such as the gravitational potential and the 
bulge mass are also explored, but these are found to be inferior when compared 
with both the bulge gravitational binding energy and bulge velocity dispersion 
predictors, for the full galaxy sample. 

Subject headings: black hole physics — galaxies: bulges — galaxies: fundamental 
parameters — galaxies: nuclei 

1. INTRODUCTION 

The presence of supermassive black holes (SMBH) (10^ — lO^^M©), in the centers of vir- 
tually all galaxy bulges, has become widely accepted in recent decades. With the acceptance 
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of the prevalence of SMBHs in the universe has come an effort to understand them: how 
they are formed, how they evolve, and how they relate to their host galaxy. The combination 
of HST spectroscopic and imaging data of the centers of nearby galaxies with ground-based 
imaging and spectroscopy at large radii, has made accurat e meas u rements of S MBH masses 
in nearby galaxies possible (e.g. iKormendy fc Gebhardtl (120011 ): iKormendyl (120041 )) for a 
large enough sample to begin to understand the SMBH-host galaxy connection. 

Given the close physical connections between the SMBH and the host galaxy bulge 
(where 'bulge' in Sections [U - El refers to either the hot, spheroidal component of a spi- 
ral/lenticular galaxy or to a full elliptical galaxy) it is possible to find host galaxy (bulge) 
characteristic parameters which may be used to predict the SMBH mass. Previous studies 
have shown that the SMBH mass (M,) is well-predicted by both the stellar bulge luminosity 
and by the associated stellar bulge mass; M, is (roughly) linearly proportio nal to both the 



bulge luminosity in the visible and NIR waveb a nds and to the bulge ma s s (IDresslei 



2003 



Kormendvl Il993l: iKormendv fc Richstone 



1995 



Haring fc Rixl |2004J ) , and I Graham et al.l (120011 ) has shown that galaxies with more 



McLure fc Dunlod l2002l : Marconi fc Hunt 



1989 



concentrated bulge light have larger SMB H masses. The SMBH mass is also well-predicte d 
by the bulge stellar vel ocity dispersion, a. jGebhardt et al. 2000b : Ferrarese fc Merritt 2000) 



such that M, oc (T^'°^ ( Tremaine et al. 20021 ). Each of these predictors {Mhuige, Lh^ige, and 



a) is physically related to each other through the fundamental plane and the virial theorem, 
and each is a physical probe of the mass distribution of the bulge. Understanding which of 
these correlations is the most fundamental would provide insight into SMBH formation and 
evolution mechanisms. 

Given the close physical relationship between these predictors [a, Mbuige, Lbuige], the 
uncertainties associated with both the SMBH masses and with the host galaxy parameters, 
and the small (< 40) number of nearby (< 110 Mpc) SMBHs with dynamically measured 
masses available in the samples used to determine these relationships, it is not unexpected 
that none of the predictors stand out as being the clear best, or most fundamental, predictor of 
SMBH mass. All of the predictors show a (statistically) similar scatter of SMBH mass about 
the predictive relationship of ~ 0.3 dex. iNovak. Faber. &: Dekell (120061 ) applied a rigorous 
statistical analysis to these galaxy property predictors from both the perspective of a Theorist 
(seeking the tightest correlation, i.e. the most fundamental predictor is that with the smallest 
residual variance) and from the perspective of an Observer (seeking the galaxy p r operty 
which can provide the best estimate of the SMBH mass). iNovak. Faber. fc Dekell (120061 ) 
determined that neither of these questions could be fully answered with the current small 
sample sizes. The probability distributions of the residual (intrinsic) variance ( Theorist) and 
the probability distribution of the uncertainties associated with the predicted SMBH masses 
( Observer) both overlapped to such a degree that it was neither possible to state which host 
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galaxy property has the most fundamental connection with the SMBH mass, nor to state 
which observational property would make the best prediction of the SMBH mass. Given the 
limited number of bright, nearby galaxies, and the demise of STIS, substantially increasing 
the sample to better address these questions, in the near future, is difficult. 

In this paper, we approach the question of what is the best predictor of the SMBH 
mass by examining two different questions: is there either a multivariate predictor or a 
previously unexplored host galaxy property which is bette r at p red icting the SMBH mass 



than those properties discussed in iNovak. Faber. fc Dekell (120061 )? iNovak. Faber. fc Dekel 



( 120061 ) found that bulge velocity dispersion, bulge mass, and bulge luminosity are all equally 
good at predicting the SMBH mass. If we assume that bulge velocity dispersion is the best 
predictor of SMBH mass, since it is not dependent on bulge-disk decompositions for non- 
elliptical galaxies, would a multivariate fit combining a and a second, or third, parameter 
be a better predictor of SMB H mass than M,{a), i.e. one with lower scatter? Second, 



Novak. Faber. fc Dekell (120061 ) comp ared bulge velocity d ispersion, bulge mass, bulge lumi 



nosity and the concentration index ( Graham et al.|[200ll ). but is there a different, physical 



host-galaxy property which is a better (or given the sample size statistically equivalent) 
predictor of SMBH mass? The motivation underpinning these questions is to gain a bet- 
ter understanding of the host galaxy-SMBH symbiosis (a Theorist's motivation), but the 
approach employed is that of an Observer; the best predictor is determined based on the 
relative predictive strengths as determined by comparing the residual scatter in the fit, not 
by comparing the intrinsic uncertainties. (However, while we approach this question from an 
Observer's framework, some of the galaxy properties explored, such as gravitational binding 
energy, are not properties which can be directly "observed"; geometrical assumptions and 
modehng are required.) For all of these calculations we use a single set of data and modeling 
assumptions which allows for direct comparisons of the relative predictive strengths. 

The paper is divided into four main sections: a discussion of the data (Section [2]), 
a description of the modeling employed to determine the host galaxies properties such as 
gravitational binding energy (Section[3]), an exploration of SMBH mass predictors (SectionH]) 
and a summary and concluding discussion (Section [5]). In this final section, we assess the 
relative merits of the predictive fits, and an argument is made for why the most fundamental 
predictor of SMBH mass, for elliptical galaxies, is the gravitational binding energy (M, oc 
Eg'^), although the bulge velocity dispersion (without additional parameters) may be more 
easily implemented, particularly in the case of spiral and lenticular galaxies. Additionally, 
there are two appendices: the first. Appendix [Al explores the selection of the host galaxy 
parameters used in Section 12.21 and the second. Appendix [Bl explores the need, or lack 
thereof, for multivariate or log-quadratic predictive relationships for the velocity dispersion, 
gravitational binding energy, gravitational potential and bulge mass predictors. 
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2. DATA 

In order to probe the connection between the host galaxy bulge and the SMBH mass, 
bulge properties (gravitational potential, gravitational binding energy, and bulge mass) are 
determined for a sample of nearby galaxies. The calculations use data taken from the pub- 
lished literature combined with geometrical assumptions and standard mathematical formu- 
lae, as describe d in S ection [31 The initial sample consists of the 30 external galaxies in 



Tremaine et al.l (120021 ). a mix of spirals, lenticulars (SOs) and ellipticals, which have black 



hole masses determined either by stellar dynamics or by gas/maser kinematics. 



2.1. Surface Brightness Profile 

The primary data required to calculate the stellar bulge mass, gravitational potential 
and gravitational binding energy are the galaxy major-axis stellar surface brightness (SB) 
profiles from which the stellar mass profiles are generated. This is implemented by calculating 
the deprojected luminosity density from the SB profile, using the methodology (and code) 



from iGebhardt et al.l (119961 ) , and then scaling the luminosity density profile by a mass-to- 
light ratio to produce the major-axis mass profile. The ideal galaxy SB profile for these 
calculations would consist of a combination of HST photometry at small radii and ground- 
based photometry at large radii, well-matched in the region of overlap, and covering the full 
extent of the major axis, in the same wavelength band as the Thpc mass-to-light ratio. For 
spiral and lenticular galaxies there is the added requirement of a well-constrained bulge: disk 
decomposition. 

The surface brightness profiles used here, given in Table [H are a combination of data 
from different literature sources, obtained using different instrumentation (both different 
telescopes and different detectors), different wavelength-bands, and different reduction tech- 
niques and corrections, along different axes; these have been combined and homogenized here 
to produce a single major-axis profile in the desired wavelength-band. In order to combine 
the data, the most compatible radial subsets of the full SB profile were selected from each 
literature source, with visually discrepant points discarded, and with no new corrections 
(K-corrections, local and/or external galaxy extinction, seeing effects and/or deconvolution) 
applied or removed. For elliptical galaxies for which both data points and a fit (Nuker-law, 
Sersic, de Vaucouleurs, etc.) were available, the data points were selected in preference to 
the fit; the fit is the result of a mathematical approximation to the real data structure, and 
although it results in a smoother SB profile, if the discontinuities are the result of physical 
structural components, they will only be imperfectly accounted for in the fit. For the spiral 
and lenticular galaxies, the fit describing the bulge light, as obtained from a literature-based 
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bulge-disk decomposition (denoted as Fd in Table [T], with the corresponding reference to the 
decomposition), is used instead of the observed combined-light profile. All of the decomposi- 
tions assumed a de Vaucouleurs profile, except for NGC 4342, which utilized a combination 
of an exponential nuclear and outer disk and a Nuker-law fit to the bulge along the minor 
axis. Galaxy light was considered to come solely from a bulge or a disk; light originating from 
bars or other features were included within these two components. An example combined 
SB profile, along with the deprojected luminosity density and its derivative (which indicates 
discontinuities in the profile), all as a function of radius, is shown for the elliptical galaxy 
NGC 6251 in Figured! this illustrates the relatively smooth connection between data from 
different origins. 

There were two corrections applied to the data: a correction to the major axis, and a 
correction to a common wavelength-band. The first correction was to rescale the SB profile 
to the major axis. Published SB profiles along both the galaxy major axis and along the 
isophotal major axis were taken to be equivalent, although this may not be strictly true in 
boxy galaxies or in galaxies with a strong isophotal twist, and no correction was applied to 
these data. For SB profiles given as a function of geometric mean radius (r = ^/ac), the 
corresponding major axis radii, for each point in the SB profile, are described as r(l — e)"'''^. 
The ellipticities were taken as a function of radius (denoted by the symbol P in Table [1]) 
when available; otherwise the ellipticity value at the nearest radial point in the profile, or an 
average value for the entire profile, was adopted. 

The second correction was a conversion of the SB profile segments into a common 
wavelength-band, the wavelength-band of the T^pc mass-to-light ratio, prior to combining 
the segments. Given that optical filters can differ from telescope-to-telescope, the published 
color terms resulting in the smoothest combined SB profiles were selected, in two cases 
(see Table [1]) requiring additional offsets to match the individual SB profile segments more 
smoothly. For direct comparisons of Jg, the V-band was selected as the default, and, for each 
galaxy, if necessary, a second (V-band) profile was generated. The color terms used are given 
in Table [TJ When available, color terms as a function of radius (denoted by P in the table) 
were used. When this information was unavailable, color terms as a function of aperture (A 
in the table) were used, or an assumption of constant color for the entire galaxy was made; 
most galaxies showed only a small (if any) gradient in color. For lenticular galaxies, if no 
separate color terms for the bulge and disk could be located, a single value was used; in 
galaxies with individual bulge and disk colors (e.g. NGC 3245), the colors in the bulge and 
disk were similar. 
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2.2. Host Galaxy Parameters 



In addition to the SB profile, the calculation of the bulge mass, gravitational potential 
and gravitational binding energy requires the distance to the galaxy (D), the galaxy bulge 
mass-to-light ratio (T), the bulge effective radius (-Re), the galaxy inclination {9), and the 
bulge apparent axis ratio (q). The SMBH mass (M,) is also required for each galaxy in 
order to generate the M,-predictive relationships in Section |H The values for each of these 
parameters are given in Table [2l Additionally, the effective bulge velocity dispersion (a) is 
required so that M,(cr) can be determined for each galaxy subsample and used for compara- 
tive purposes in Sections S] and El the velocity dispersion values are enumerated in Table [31 
along with the cor nputed host galaxy properties. The distance for each of the galaxies is 
taken direct l y from iTremaine et al.l (|2002| ). iTremaine et al.l (120021 ) used SBF distances from 
Tonry et al.l (120011 ) . when available; otherwise the distance was determined from the reces- 
sion velocity assuming a Hubble constant of 80 kms~^ Mpc'^ . All parameters, including the 
SMBH mass and mass-to-light ratios, are scaled to the distance (and, implicitly, use this 
Hubble constant). The SMBH masses are taken from the individual black hole modeling 
papers, along with accom panying estimat e s of t he uncertainties. Generally the distance- 
rescaled value agrees with ITremaine et al.l ( 20021 ): in the few galaxies where this is not the 
case, a notation is made in Table [2l The effective bulge stellar velo city dispersions ( a, defi ned 
as the rms dispersion within a 2Re slit aperture) are taken from Tremaine et al.l (120021 ) for 
all galaxies except NGC 4258, which was taken from lSiopis et al. J2006h . 



The dynamical mass-to-light ratio (T) is used to convert the stellar luminosity density 
into a stellar mass density. The strong dependence of the computed mass, gravitational 
potential and gravitational binding energy on the mass-to-light ratio makes it crucial to have 
the most accurate T-value possible. The implications of selecting an "incorrect" value, and 
the range of values quoted in the literature, are explored in Appendix lAl The best mass 



to- ligh t ratio values, available in the literature, are those determined by ICappellari et al. 



(120061 ) using extinction-corrected NIR photometry and integral-field-spectroscopy. These are 
an improvement on other mass-to-light ratios both because of the improved spectroscopy, 
and because observations at redder wavebands minimize the impact of dust; the mass-to- 
light ratios determined at blue wavebands are probing the star formation history as well 
as the mass. However, these superior mass-to-light ratios are only available for one-third 
of the total galaxy sample; therefore, for this subset of the sample galaxies, two mass-to- 
light ratio values are enumerated. For every s ample galaxy, th e best mass-to-light ratio 
available in the literature, excluding the value in ICappellari et al.l (120061 ). is determined and 
denoted as the best-pre-Cappellari value, Tbpc (Table [2l column 6). For galaxies in which 
T was determined in the process of measuring M, (using well-determined spectroscopy and 
photometry from both HST and ground-based observatories), this value was adopted as 
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T 



hpC- 



For galaxies in which T was not provided along with M,, Tbpc was taken from 



Tremaine et al.l (|2002| ). if available, and otherwise from another literature source (Table [2) 
column 7). The wavelength band of T^pc is taken as the default for the galaxy; the SB 
profiles were converted to this wavelength-band prior to calcu lating the luminosity d ensity 
profile. If the mass-to-light ratio for the galaxy was provided in ICappellari et al.l (120061 ) then 
this T, rescaled to the Thpc-default- wavelength band, is given as Tcap in Table [2] (column 
8). The best mass-to-light ratio, T^est = ^, used in the galaxy modehng, is taken to be Tcap 
if it exists; otherwise, TbpC is utilized. As discussed further in Section HJ mass-dependent 
predictors of M, have lower residual scatter when using the Tcap-values than when using 
the Tfcpc-values, for galaxies in which both are available. It is assumed throughout that the 
mass-to-light ratio is constant with radius (T(i?) = constant). This assumption was also 
made by investigators when determining M,; authors who investigated a radial variation in 
selected sample gala xies generally did not fin d a large variation in mass-to-light ratio. For 



example, in IC 1459 ICappel 



ari et al.l (2002 ) estimate a decrease in T of 15% per decade 



in radius, and in NGC 3379, iGebhardt et al.l (j2000al ) estimate T may be 25% higher in the 
innermost regions than in the outermost regions of the galaxy. 

The observed bulge axis ratio (q) and bulge inclination (6), in combination with sev- 
eral simplifying assumptions (Section [3]), are used to geometrically project the major-axis 
mass profile (Section 12. ip into a 3-dimensional mass profile, assuming the bulge is an oblate 
spheroid. The inclination and axis ratio are also used to correct the observed SB profile 
and luminosity (and mass) density profiles to an edge-on orientation using a multiplicative 
rati o of the observed (q) and intrinsi c , edge -on, (p) axis ratios; [(psin6')^ = ~ {cosOY 
e.g. iRichstonel Jl979h : iGebhardt et al.l (119961 ) . The predictive fitting functions are based on 
values determined from the edge-on galaxy orientation, although the results are qualitatively 
unchanged if this correction is omitted. As in the case of the mass-to-light ratio, there is 
often a range of q- and ^-values present in the literature; the range of these values and the 
implications of selecting an "incorrect" value for q or ^ are discussed in Appendix \^ 

The observed axis ratios, culled from the literature, are assumed to be constant as a 
function of radius for all galaxy bulge radii. If multiple values for the axis ratio were quoted 
in the literature, the most prevalent literature-value, or an average of the literature values, 
as given in Table[2J columns 10-11, was selected as the best value. For the elliptical galaxies, 
the best axis ratios are taken from the literature either in the form of a single quoted value 
representative of the galaxy as a whole, or derived from the radially-dependent ellipticities 
associated with the SB profiles. For the spiral/lenticular bulges the published literature 
single-values are not used; such single-values are based on a combination of both the bulge 
and the disk light profile, and so are flatter than for the bulge-component alone. Instead, 
the bulge axis ratios are based on the SB radial profile ellipticities at representative bulge 



- 8 - 



radii. 



The bulge inclinations, for all galaxy morphologies, are assumed to be constant as a 
function of radius (i.e. no warping) and, in the case of spiral/lenticular galaxies, to have the 
same inclination as the galaxy disks. The inclinations used in the determination of the M, 
values are adopted, when given. If such a value was unavailable, an inclination was taken 
from the literature, and if this, too, was unavailable, the inclination was determined from 
the observed axis ratio. Inclinations are generally better-determined for the spiral/lenticular 
galaxies, for which the inclinations can be measured geometrically from the observed disk 
axis ratio, than for the disk-less elliptical galaxies. For several of the galaxies, the M,- 
determination modeling implies, or assumes, an edge-on inclination for the galaxy, while 
independ ent structura l anal yses of the galaxies imply a less-than edge-on inclination. For 
example, iBower et al.l (120011 ) find that an inclination of 90° results in the best-fitting three- 



i ntegr al model, used to determine the SMBH mass, for NGC 1023, while iDebattista et al. 



(120021 ) fin d, in an examination o f the bar of NGC 1023, that the best inclination is 66°.4±1°.2. 
Likewise, iBusarello et al.l (119961 ) claim a less-than edge-on inclination for NGC 3384, follow- 
ing their structural analysis of that galaxy. These inclination discrepancies, however, do not 
affect our conclusions, as illustrated in Appendix 1X1 

The effective (half-light) radius of the bulge (-Re) is not directly used in the calculations; 
it is used as a fiducial reference point for the gravitational potential and other radially- 
dependent functions. As with the previous bulge parameters, there were often multiple values 
for i?e in the literature. However, an incorrect value will not affect the intrinsic mass profile or 
the computed gravitational binding energy; it will only affect predictors which are evaluated 
at nRei and so have a minimal impact on the conclusions. For the spiral/lenticular bulges 
the bulge effective radius was taken from the bulge-disk decomposition literature source; the 
exception is NGC 4258 which wa s obtained direct l y from C. Siopis (private communication). 
For the 10 elliptical galaxies in iGebhardt et al.l (120031 ). the value was taken directly from 
that paper. For the remaining elliptical galaxies, the value from the literature, or from an 
unweighted r^/^ law fit to the SB profile, which resulted in the least scatter in a fundamental 
plane relationship (log Re = a log le + b log a + c) fit to the full sample of eUiptical galaxies 
was chosen as the best. For most galaxies it was assumed that the effective radius is not a 
function of color, i.e. that the color is constant as a function of radius. 
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3. Determination of the Stellar Mass, Gravitational Potential, and Binding 

Energy 

In order to determine the bulge mass, gravitational potential and gravitational binding 
energy, the input data (see Section [2]) are combined with well-constrained mathematical 
functions and several simplifying assumptions about the host galaxies. The first assumption 
is that the galaxies are smooth and featureless; however, while this may be the case for 
NGC 221, it is certainly not the case for the rest of the galaxies which include, among other 
features, nuclear star clusters, AGNs and cores, ionized gas disks and shells, dust arranged 
in clumps, filaments, rings and disks, outer and inner disks which are sometimes warped, 
and bars and jets. (For details on the individual galaxies, see the BH modeling and SB 
papers previously referenced.) A second assumption is that not only is the galaxy uniform, 
but that it is an axisymmetric oblate spheroid with constant mass density along concentric, 
isodensity, oblate spheroidal shells with no radial variation in ellipticity; if galaxies are 
triaxial or have isophotal/isodensity twists with radial position, this will be unaccounted 
for in these calculations. Third, it is assumed that there is no contribution from inner or 
outer disks, bars or the dark halo. The binding energy (as well as the mass and potential) is 
considered to come solely from the mass associated with the visible stellar light; the binding 
energy could be much greater if the galaxy is immersed in a massive dark halo. While 
ignoring the dark halo, which can only be inferred, rather than directly measured, will have 
an impact on all galaxy types, ignoring the presence of visible disks in the spiral and lenticular 
galaxies may have a substantial effect on these galaxies, weakening the predictive power of 
the calculated parameters relative to an elliptical galaxy sample. 

The computation of the bulge mass, gravitational potential, and gravitational bind- 
ing energy is implemented using well-characterized mathematical expressions under the as- 
sumption that the bulge is a smooth, featureless, oblate spheroid. For an oblate spheroid 
{a = b > c) with axes a, b and c, the radii along the major axis are projected to any ar- 
bitrary radii using the relationship that r(z/) = r^naji^ + /cqZ/^)^^-^ where = a/ (a/c)^ — 1 
and u = cos 6. (In this section, r, 6, and refer to standard spherical geometry coordinates.) 
This allows for the projection of the major-axis mass density profile (see Section [2?T|) to any 
arbitrary axis: p(r, z/) = p(r(z/),0). With this relation, the mass enclosed within a radius r 
can be calculated as 

M = p{r, u)r^drdud(f). (1) 



The ca lculation of the gravitational potential, $, is derived from equation 2-122 of lBinney fc Tremaine 



( ll987l ). under the assumption that the density does not depend on 0, to be 

$(r, = -27rG'^P/(cos z/) [r^^'+^M^ + r^Bi] (2) 
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where ^ 

Ai = / aW Pi{v)p{a,v)dvo?da 
Jo J -I 

poo pi 



/oo pi 
a"('+^) J Pi{u)p{a,u)dua^da. 



U sing; this relation for poten tial, the binding energy can be calculated, as in equation 2-19 



of iBinney fc Tremaind (119871 ) as 



E = - / / p{r,u)^{r,i')r'^drdud(l). (3) 

2 J(j) Ju Jr 

Ou r computational c ode utilizes multivariable, Gaussian quadrature integration tech- 



niques (jPress et al.l (119921 ) routine qgaus) combined with simplifying assumptions about the 
radial extent and structure of the mass-density profiles. First, the mass profile (Section [2.11) 
is specified at discretely sampled points between arbitrary limits. It is assumed that the 
profile is intrinsically smooth, and the space bet ween points is interpolated using a cubic 



spline algorithm: spline from iPress et al.l (119921 ). Second, the mathematical expressions 



(equations [D - [3]) integrate over the spatial variable (r) from zero to infinity. It is assumed 
that the galaxy mass physically extends only from an innermost radius {Rmm 7^ 0) to an 
outer limit {Rmax 7^ cxd), and outside of these limits the mass density is taken to be zero. 
Third, the value of Rmin from the observationally-based mass profile reflects the resolution 
limits of the observation, not the inner cutoff of the galaxy mass. The mass profile was 
extrapolated inward to 3 Schwarzschild radii, using an r — logp unweighted quadratic fit 
to the innermost 10 points of the mass profile using a least squares fitting algorithm {Isqfit, 
discussed in Section H]). The computed energy for fitting-sample galaxies was unchanged if 
the extrapolation is repeated using an assumption of constant mass or if the innermost limit 
is varied to other physically reasonable values such as 0.01 pc. Fourth, the value of Rmax 
likewise reflects the observational limitations rather than the physical galaxy edge. The 
mass profile was extrapolated outward (to an extreme limit of 50-Re) using a logr — logp 
unweighted quadratic fit to the outermost 10 points of the mass profile using Isqfit. The 
computed energy for fitting-sample galaxies was unchanged both if the outermost limit is 
varied to any value beyond a few Re and if it is assumed that the observational cutoff to the 
mass profiles corresponds to the physical edge of the mass profile. Finally, the integration 
only goes out to a Legendre polynomial of order 4. Tests using models with homogeneous 
spheres, Plummer density profiles, and Satoh density profiles indicate that these selections 
provide sufficient spatial coverage and that the expected potentials are produced. The final 
calculated values for the bulge stellar mass enclosed by radii R, the stellar gravitational 
potential evaluated at radii R, and the stellar-based gravitational binding energy, for all 30 
galaxies in the original sample, are given in Table [31 
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PREDICTORS OF BLACK HOLE MASS 



Utilizing the host galaxy parameters described above, predictors of the SMBH mass, 
primarily in the form of power laws, were determined using a variety of algorithms. The 
host galaxy-SMBH mass predictors are of the form 



F = ^ aiXi + d 



(4) 



where Xi are the host galaxy parameters, Y = logM,, and Xj = logXj. The host galaxy 
parameters were each normalized by a value near the mean (in log-base- 10 space) of the 
30 galaxy sample (Table [3]) • This norm alization remove s the c ovariance between the zero- 



point (d) and the slope, as discussed in iTremaine et al 



rithm is Isqfit (based on formulae/methods in iBevingtonI (119691 )): this program calculates 



20021) ■ The primary fitting algo- 



the weighted-least-squares-minimized fit to a multivariable equation of the form of equa- 
tion HI weighting points only by the inverse-square of the uncertainties in Y associated with 
each data point, not by the uncertainties in Xj. In addition to the coefficients (oj, d), the 
algorithm provides coefficient uncertainties {6ai, 6d) using standard least-squares-fitting for- 
mulae, e.g. iPress et al.l (119921 ). The reduced chi-squared for this function is determined such 
that 

1 



DOF 



N 



-2 



[y;-(5^a,x,, + rf)fj 



(5) 



where DOF is the number of degrees of freedom a nd ej = efy (e,, is the total uncertainty 
in logM,). The second fitting algorithm is fitexy (jPress et al.lll992l ): this only allows for 



the fitting, via least-squares-minimization, of one parameter {i = 1, equation Hj), but, unlike 
Isqfit, it includes uncertainties in both Y and Xi. The reduced chi-squared for this function 



is determined by equatio n [5] where j = 1 and 
dependent variable Xi 



[ex is the uncertainty in the in- 
Novak. Faber. fc Dekell (l2006l ) found this algorithm for determining 
SMBH-predictive relationships to be the best, i.e. the most e fficient and least b iased among 



a set of seven algorithms explored. The third routine, medfit (jPress et al.lll992l ). henceforth 
robust, also fits to only one parameter, but utilizes absolute-deviation-minimization, and 
does not use the uncertainties in either Y or Xi. 

The measured uncertainties on M, are not symmetric; therefore, the calculated fits are 
dependent on the method of symmetrization employed. The default procedure (referred to 



yOhigh 



+ e 



yOlon 



where e^o is the 



as Avgerr) is to average the uncertainties such that eyo = 0.5 [e 
observationally-based uncertainty. The second method (Recent) involves recentering M, be- 
tween the upper and lower observational limits, such that Y = 0.5[{Y + €yo,^-g^) + (F — ^yoi^^)]- 
For all fits, the Avgerr and Recent fits are equivalent, unless specifically noted in the text or 
tables. 



- 12 - 



Additionally, although the fits minimize x^? even with the tre atment of error i n bot h 



variables, the resulting Xr exceeds 1.0 (see equation [5]). Following iTremaine et al.l (120021 ) 



we assume that there is an additional variance (eyj„) due to cosmic scatter and errors in 
the independent variable, and account for it by adjusting the error in Y = logM, to = 
^yo + ^yin- The intrinsic uncertainty is selected in order to always obtain = l-O- The 
predictive fit is then calculated using three weighting methods. In the first (OBS), it is 
assumed that there is no intrinsic error (e^ = eyo); thus Xr l-O- the second method 
(INT), the observationally measured uncertainties are ignored (e^ = eym)- In the final, 
default, method {OBS+INT), it is assumed that the total uncertainty is a combination of 
the observed and intrinsic uncertainties, as given above. For all fits, the three methods are 
equivalent, unless specifically noted in the text or tables. Throughout, the term weightcent 
will be used to refer to the selection of the combination of Y-centering [Avgerr or Recent) 
and Y- weighting {OBS, INT or OBS+INT). 

The final weighting option is in the adopted uncertainty in Xi (e^,) for the fitexy algo 



rithm. For the bulge velocity dispersion this is set to the value adopted in iTremaine et al. 



(120021 ): ex = 0.021 (5%). For the gravitational binding energy, gravitational potential, and 
bulge mass, values of = 0.1 (26%) and = 0.3 (100%) are explored to see if the fits are 
substantially different from the default case, = 0.0. 

Throughout the paper the terms strength, stability, and best are used to describe the 
fits. The strength refers to the amount of scatter {cfy) about the predictive relationship, 
where 



^fy=\> ^Y^-0 . a,X,, + d)] 1 DOF-'. (6) 



Less scatter implies a stronger fit or predictor. A stable slope/fit is defined to be one for 
which the variations in the slope/fit, with changes in weightcent, fitting algorithm and/or 
morphological selection, are statistically insignificant. The strongest and the most stable 
predictor is referred to as being the best predictor. 

The stability of the fits i s further examined for each sample using a bootstrap calculation. 



A bootstrap calculation, e.g. iPress et al.l (119921 ). randomly draws N points, with replacement 



from an initial set of N data points, resulting in a 'new' data sample with 63% unique data 
points. This method indicates whether the fit based on the original sample is being unduly 
infiuenced by the presence of specific galaxies, and provides insight into the underlying 
distribution of the best-fit coefficients. These tests are only intended as a consistency check 
on the slopes and uncertainties based on the full, original sample of N data points, and 
to further illustrate the robustness of the results. For each galaxy sample, the bootstrap 
is run for 1000 'new' samples using both the default weightcent, Isqfit fitting algorithm 
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and the robust algorithm. The value of for the Isqfit fitting is selected to produce 
xl = 1.0, assuming that only unique data points in the 'new' sample count towards the DOF 
(counting duplicated galaxies in the DOF as individual galaxies had a minimal impact on the 
resulting slope except for the 8B sample.) The mean, median, average deviation, standard 
deviation, skew and kurt osis of the di s tribut ion of values are determined using the selip and 
moment algorithms from iPress et al.l (119921 ). for the slope (a), zero-point (d), scatter {(Jfy) 



and required intrinsic uncertainty in F = logM, (eyin). 



The scatter (cfy) and the Snedecor F-test (ISnedecor fc Cochranlll989l : iDixon fc Massey 



19691 ) are used to statistically assess the relative goodness- of- fit for the predictive relation- 



ships (for a given sample) . Usin g as a comparison is not f easible since = 1 for the least 



square fitting. The F-test ratio (ISnedecor fc CochranI Il989l : iDixon fc Masseyl Il969l ) , 



(7) 



[where F^^ > 1.0 implies a better fit than M,((t)], compares the square of the residual 
scatter of points around each of the two predictive relationships and determines whether the 
difference in scatter is significant given the number of degrees of freedom (DOF) in each of 



the re lationships. This determination is made using the published tables in lDixon fc Massey 



(119691 ) and results are abbreviated as follows: NSS-the fits are not different at the 75% 
level (approximately l-a), 75SS - the fits are different at the 75%-90% level, 90SS - the fits 
are different at the 90%-95% level, 95SS - the fits are different at the 95%-97.5% level and 
99SS - the fits are different at greater than the 99% level. (For ease in using published 
tables, assumptions that 13 D0F~12 DOF {15E) and 21 DOF~20 D0F~19 DOF {23gal) 
were made; these assumptions do not change the stated results.) Throughout, a mention 
of 'minimal improvement' implies that the scatter is lower than M,((t), but not at the 75% 
probability level. The abbreviations and symbols defined in this section are summarized in 
Table H 

The remainder of this section will be structured as follows. First, the selection of 
galaxies used in the four fitting samples will be explained, along with a caveat related to 
using the spiral/lenticular galaxies for fitting (Section 14. ip . This will be followed by the fits 
for the bulge velocity dispersion predictor (Section 14.21) . the multivariate fits constructed 
from a combination of a, le and Re (Section 14.31) . the gravitational binding energy predictor 
(Section|13D, the gravitational potential predictor (Section 14. 5p and the bulge mass predictor 
(Section l4.6p . Each of the single- variable predictors (Sections 14.21 14.41 - 14.61) . will begin with 
a discussion of results from the literature (if available) and the motivation/methodology for 
fitting this galaxy property. This will be followed by the fits for each of the samples and 
any notable fitting caveats, along with a discussion of the stability of the fit in terms of 
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the galaxy morphology, the fitting algorithm, and the weightcent selection, and in terms of 
the bootstrap results. Finally, plots of the fits and their residuals will be examined, and 
comments on the necessity, or lack thereof, for either a multivariate fit or a log-quadratic fit 
will be made. Additionally, for all but the M,{a) predictor, the strength of the fit compared 
with the M,((t) fit is assessed by looking at the F-test ratios. 



4.1. Sample Selection 



The initial sample of 30 external galaxies (ITremaine et al.l 120021 ) contains 17 elliptical 
galaxies, 10 lenticular galaxies and 3 spiral galaxies; of which 2 elliptical, 2 spiral and 3 
lenticular galaxies are excluded from future discussion and fitting for the following reasons. 
The elliptical galaxy NGC 2778 is rej ected because it is an outher in most relations, has 
been noted to contain a stellar disk (IGraham et al.l 1200 ll : iRix. Carollo. fc Freeman! Il999l ) 
like a lenticular, and has a relatively large observational uncertainty (although its inclusion 
would not change the final fits because of its associated low relative weight). The elliptical 
galaxy NGC 4564 and the lenticular galaxy NGC 3384 are both rejected because of their 
relatively small observational uncertainties (high weight) combined with the fact that they 
are marginal outliers in some of the relations (particularly those involving for NGC 4564); 
their inclusion results in noticeably different fits. The spiral galaxies NGC 1068 and NGC 
4258 are rejected because they are complicated spiral galaxies containing AGNs; obtaining a 
correct bulge- hght-only surface brightness profile, equivalent to elliptical galaxies, is difficult, 
and these galaxies are usually outliers in calculated fits derived from the SB profiles. The 
lenticular galaxy NGC 4342 is rejected because it, too, has a complicated surface brightness 
profile resulting from the combination of a prominent inner and outer disk; furthermore, this 
gal axy has been n o ted to be an outlier in other relationships, such as the bulge mass predictor 
by iHaring fc Rixl J2004h. The lenti c ular g alaxy NGC 4742 is rejected because the SMBH 
mass, tabulated in ITremaine et al.l (120021 ). has not bee n described in a separ ate, detailed 
analysis paper. The Milky Way, the remaining galaxy in ITremaine et al.l (120021 ). is not used, 
despite its precise SMBH mass measurement, because there is no way to reproduce the 
surface-brightness-pro file-driven calculations in a manner identical to the external galaxies. 
Tremaine et"aD J^OoJ identifies 9 of the galaxies (the Milky Way, M31, NGC 1068, NGC 
2778, NGC 3115, NGC 3379, NGC 5845, NGC 4459 and NGC 6251) in the full-31 galaxy 
sample as being 'questionable'; however, not all of these are excluded from this analysis, based 
on our adopted selection criteria. The use of even more stringent selection criteria would 
render the samples too small to make any st atistically mean i ngful comments. However, 
galaxies rejected from inclusion in the initial ITremaine et al.l (120021 ) 31-galaxy sample as 
having unreliable M, are excluded from the fitting here. 
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The final-fitting-sample consists of 23 elliptical, spiral and lenticular bulges {23gal sam- 
ple); this full sample is then subdivided into 3 final-fitting-subsamples based on morphology 
and availability of data. The full, 23gal sample is first subdivided into a sample contain- 
ing the 15 pure ellipticals {15E) and a sample containing the 8 remaining spiral/lenticular 
bulges {8B). Finally, there is a subsample of the 15E galaxies cons i sting of 8 elliptical galax- 



ies {Cap8) for which mass-to-light ratios from ICappellari et al.l (120061 ) are available (see 



Section [2l2l) . These samples are summarized in Table |H 

A caveat related to fitting with the spiral/lenticular galaxies (present in both the 8B and 
23gal samples) will be noted before explicitly discussing individual predictors in the following 
sections. The SB sample exhibits problems when weighting (in Y) only by the observational 
uncertainty {OBS) which are not found for the elliptical sample. These problems likely 
stem from either incorrect eyo values for some of the galaxies (which results in incorrect 
weights) or from incorrect host galaxy parameters. The calculated OBS-weighted fits for 
the spiral/lenticular galaxies have large and scatter and have significantly different slopes 
than when /A^^T- weighting is used. For example, the OiJ^*- weighted M,{a) relationship for 
8B: log (M./Mq) = (1.82 ± 0.88) log {a / {200kms-^)) + (7.81 ± 0.11) {scatter = 0.53) bears 
no resemblance to the accepted slope (which is near 4.0) and is an extremely poor fit for the 
elliptical galaxies. This problem is prevalent for all of the predictive relationships examined; 
adding the excluded 5 bulges does not solve the problem, and removing additional bulges 
results in too small a sample to make any statistically meaningful comments. Although the 
OBS weighting is not being used for the final fits, the observational uncertainties do factor 
into the OBS+INT final fits, and so the predictors based on both the 8B and the 23gal 
samples of galaxies may be less reliable than fits based solely on the elliptical galaxies. 



4.2. Bulge Velocity Dispersion 

The first M, predictor examined is the bulge stellar velocity dispersion (a); this rela- 
tionship is used as a point of comparison for the predictors in the future sections. Table 
presents the predictive fits (equation Hj) as follows: for each sample {15E, 23gal, SB, CapS) 
and fitting algorithm used {Isqfit, fitexy, robust), using the previously specified default weight- 
cent options {Avgerr, OBS+INT), the slope (a) and the zero-point (d) of the fit, along with 
the intrinsic uncertainty in y = logM, (e^m) required to produce Xr = 1; the uncertainty 
in X = logo" (e^) applied, the maximum (amax) and minimum (a^m) slope obtained with 
any of the viable weightcent combinations, and the scatter (equation [6]) are given. The 
quantities/abbreviations are as previously specified (see Table H]). 

The fits for the four samples are, as one would expect, statistically equivalent to the 
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fit from tfie iTremaine et al.l (120021 ) 31-galaxy sample; tlie four samples are subsets of the 
31-galaxy sar aple and the vel o city d is persion values f o r all fi tting-sample galaxies are taken 
directly from iTremaine et al. J2002h . ITremaine et al.l (120021 ) finds a best-fit relationship for 
the 31-galaxy sample of log (M./Mq) = (4.02±0.32) log(a/(200fcms-i)) + (8.13±0.06), with 
an intrinsic uncertainty of 0.27 (or 0.23 ± 0.05 dex using a maximum likelihood estimate). 
The slopes (Table [5j column 3) for the four (sub)samples are all compatible with this result, 
within the stated uncertainties, as are the zero-points and requisite intrinsic uncertainties in 
Y ( £,„:„ ,) . It may be noted that the fitted slopes are con sistently lower t 



( 120021 ) slope of 4.02 ± 0.32. This is not unexpected; ITremaine et al 



lan the Tremaine et al. 



(l2002h states th at the 



4.02 slope "may slightly overestimate the true slope by 0.1-0.3" . ITremaine et al.l ((2002) made 
this comment after examining 4 Tremaine- subsamples which (1) excluded the Milky Way 
(a = 3.88 ± 0.32), (2) excluded high [a > 250kms-^) d ispers i on galaxies (a = 3 -77 ± 0.49), 
(3) excluded all but the 10 galaxies in iPinkney et al.l (120031 ): iGebhardt et al.l (120031 ) (a = 
3.67 ± 0.70) and (4) removed 9 'questionable' galaxies: the Milky Way, M31, NGC 1068, 
NGC 2778, NGC 3115, NGC 3379, NGC 5845, NGC 4459 and NGC 6251 (a = 3.79 ±0.32). 

The morphological division of the full {23gal) sample into the elliptical {15E, Cap8) 
and spiral/lenticular bulge {8B) subsamples illustrates a difference: the elliptical samples 
consistently have a larger slope than the bulge sample. This implies that altering the ratio of 
the spiral/lenticular :elliptical galaxies in a combined sample of galaxies will alter the slope 
of the combined sample. It should also be noted, however, that while the bulge {8B) slopes 
are lower, they are consistent with the elliptical slopes within the (large) stated uncertainties 
of the slope. 

The morphological division of the 23gal sample also illustrates that the (statistically) 
strongest fitting sample is the full-elliptical {15E) sample (also see Section 14. ip . This is 
evident both in the resultant scatter from the fits {cTfy) which is lowest for the 15E sample, 
even when compared with the full-galaxy sample {23gal), and in the required intrinsic un- 
certainty in log M, [eyin) which is consistently lower than for any other sample. The CapS is 
(statistically) weaker because it has only about half of the number of points. The 8B (bulge) 
sample is weaker both because of the small number of points, and because it is intrinsically 
weaker (as illustrated in comparing the scatter with the same-sized Cap8 sample). This is 
likely due to complications related to the increased amount of galaxy structure. Finally, the 
23gal full-galaxy sample is (statistically) weaker because of the presence of the 8B galaxies 
which are either fundamentally different (as illustrated by their lower slopes) and/or appear 
different because of imperfections in the data. 

The replaceme nt of the SMBH mass for NGC 224 (4.5xlO^M0) with the higher value of 
1.4xlO^M0 given in Bender et al. J2005h (see Table [2]), and the replacement of the velocity 
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dispersion for NGC 2787 (140 kms"^) with the higher value of 200 kms~^ suggested by 
hterature measurements of the central velocity dispersion (see Table [3]), does not change the 
conclusions drawn from the fits, as shown in Table O The slope of the 23gal least-squared 
[Isqfit) fit is essentially unchanged, while the slope of the 23gal robust fit is substantially 
lowered (3.49 compared with 3.70), but is consistent with the originally stated slope-range. 
The slope of the 8B Isqfit fit is even lower than originally stated, 3.10 ± 0.96 compared 
with 3.27 ± 0.77, while the slope of the 8B robust fit is higher, 3.66 compared with 3.34, 
but still falls below the elliptical galaxy sample slope. This illustrates both the relative 
instability of the fitted slopes for the 8B sample, and the dichotomy between the elliptical 
and spiral/lenticular bulge sample slopes. For all of these fits, the revised parameters result 
in larger values for the scatter, and in larger requisite uncertainties in logM, {eym) for the 
least-squared fitting, and so the revised parameters are not generally preferred over the 
original values. 

In addition to providing information about the predictive fits based on morphological 
divisions. Table [5] also provides information about the importance of the fitting algorithm. 
First, it is evident that, for a given sample, including (or excluding) the uncertainty in 
X = logo" {e^ = 0.021) has little impact on the resulting fits. The slope increases by only 
0.04 (1 %) when it is included, well within the stated uncertainties, and eym, the requisite 
intrinsic error in Y = logM,, decreases by only 0.01. Second, the slopes and zero-points 
resulting from the least-squares-fitting algorithms are consistent with the results from the 
robust fitting algorithm; there is no change in the conclusions or in the stated morphological 
trends. 

Altering the weightcent selections, illustrated in the stated maximum {amax) and min- 
imum (a^m) slopes, does not change the stated conclusions related to morphological type; 
elliptical galaxy samples persistently exhibit higher slopes than the spiral/lenticular bulge 
galaxy sample. The weightcent selections can result in slope variations, relative to the de- 
fault, of up to 11.5% {fitexy, SB), but these variations are well-within the stated default-slope 
uncertainties. The examination by morphological type also indicates that the default slope 
for ellipticals is at the high end of the range suggested by a^ax and a^jn, while the bulges 
are at the low end, maximizing the apparent discrepancy in the slopes. However, while the 
full ranges bring their slopes closer together, the elliptical galaxy sample slopes are still, on 
average, higher, and higher for every specific weightcent selection. 

The stability of the fits, indicated by the statistically negligible variation in slope for each 
sample with changes in fitting algorithm and weightcent, is reinforced through the bootstrap 
calculations (Table [6]). The value of is the number of 'new' data sets (out of 1000) for 
which an eyin resulting in = 1.0 could be determined; samples in which there were too 
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few (unique) points for this to occur, even when eyin = 0, were excluded. The bootstrap 
exhibits large uncertainties (both Isqfit and robust) associated with both the Cap8 and 8B 
sample slopes and large values for the skew and kurtosis; this underscores the problems 
inherent in making statistical analyses based on small samples. The 8-galaxy samples are 
reduced to (typically) 5 unique galaxies here, and with this number of galaxies the slopes 
are very dependent on individual galaxy effects (and selection) rather than on a general, 
broad description of the morphological group. Given this caveat, it is notable that the bulge 
{8B) slopes are much higher in the bootstrap results; higher than even the elliptical and 
combined samples. This dramatic change is not seen in the similarly small CapS sample, 
further illustrating that the bulge-only sample is weaker when making statistical claims, 
likely owing to the strong influence of several (discrepant) galaxies on the slope. Removing 
these galaxies would result in further problems from an even more reduced sample size, 
however. The results for the larger samples, 15E and 23gal, are generally consistent with 
the original results, both for the robust and for the least-squares fitting results. 

When the best fits, using the Isqfit algorithm, for each sample (Table [5]) are plotted, 
further information about the predictive relationships can be gleaned. Looking at the bot- 
tom half of Figured (M, versus a) for ellipticals only (left), spiral/lenticular bulges only 
(center) and the combined sample (right), the bulges are moderately fit by the elliptical-only 
relationship (dotted), while the ellipticals are not well-fit by the (lower) bulge slope (dashed 
line). Figure [3] (ellipticals in top panel, bulges in middle panel and combined sample in 
bottom panel) plots the (Isqfit) fit residuals against the bulge velocity dispersion (left), the 
bulge effective radius (center), and the intensity at the bulge effective radius (right). Neither 
in Figure [2] nor in Figure [3] does there appear to be strong evidence for the necessity of a 
quadratic, as opposed to a linear, parameterization; this is further explored in Appendix [B] 
where formal fits of a quadratic are shown to provide no improvement in the quality of the 
fit (as indicated by CTfy). There also appears to be no evidence of a pattern in the residuals 
with a, Re, or for either the combined {23gal) or for the elliptical {15E) sample. (There 
is a minimal ^S-fit-residual correlation with a when looking at only the elliptical galaxies, 
indicative of the previously mentioned poor lower slope fit for the elliptical sample). The 
bulge sample (55-filled circles), however, shows a clear correlation between the fit-residuals 
and Rf. for all three {15E, 8B, 23gal) of the predictive fits. 

4.3. The Combination of Bulge Velocity Dispersion, Radius and/or Intensity 

The M,(cr) predictor is successful for all galaxy morphologies, but the residual scat- 
ter prompts the question of whether a multivariate fit, or an alternative single- variable fit. 
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would reduce this scatter. iMarconi fc HuntI (120031 ) claim that there is "a weak, but sig- 



nificant, correlation of [M, — cTg] residuals with R^" and that a "combination of both ae 
and Re is necessary to drive the correlations between M, and other bulge properties": the 
so-called 'fundamental plane of SMBHs'. The M, — M,(cr) residuals for the elliptical {15E) 
and combined {23gal) samples, illustrated in Figure [3l do not immediately indicate such a 
combination. For the bulge-only {8B) sample, howe ver, there is good evide nce of a corre- 



lation between the residuals and R^- Examining the iMarconi fc HuntI (120031 ) data, it is the 
lenticular /spiral galaxies which are generating the appearance of a significant correlation. 
If only their spiral/lenticular galaxies are considered, there is a relatively strong correlation 
present, but looking at only their elliptical galaxies there is, at best, a minimal correlation 
between the residuals and Re- Thus, cr is a weaker predictor of M, for spiral/lenticular 
galaxies than for elliptical galaxies. 

Multivariate fits combining a, le and Re are examined, using the Isqfit algorithm (default 
weightcent) , to determine if there is an obvious, physically- motivated, predictor of M,, other 
than M,(cr). The multivariate fits (equation HI where xi = le/ho, X2 = Re/Reo, and 
xs = a/ao) for the combined {23gal), elliptical-only {15E) and spiral/lenticular bulge-only 
{8B) samples are presented in Table [71 In addition to the best-fit coefficients (a^, d), the 
intrinsic uncertainty in logM, (cyj^) required to obtain Xr = 1-0) ^^e scatter (equation [6]), 
the F-test ratio {F„y) relative to the M,{a) fit (equation [7]), and the significance (Sig) of the 
difference relative to the M,{a) fit are given. The fits, excluding M,(/e) and M,{Re), are 
illustrated in Figure S] wherein the SMBH mass predicted by the fit (based on the top panel: 
15E sample, central panel: 8B sample, and bottom panel: 23gal sample) is plotted against 
the measured SMBH mass for, from left to right, the M,((t), M,{Ie, Re,cr), M,{Ie, Re), 
M,{Re,a), and M,{Ie,cr) relationships. 

The addition of Jg or Re as a second parameter [M,(Je, cr) or M,{Re, a)] does not improve 
the quality of the fit (based on the scatter) for the elliptical or combined samples, despite 
the additional free parameter. Furthermore, the coefficient of the additional parameter (oi 
and 02 respectively) is consistent with zero for the elliptical sample and is minimal (in 
comparison with the cr-coefficient) for the combined sample. The bulge {8B) sample, which 
was indicated to require a M,{Re, a) fit by the Figure [3]residuals, does result in predictive fits 
that are significantly better than the M,(a) predictor: M, oc /-i-02^4.68 ^ ^o.94^2.90_ 

However, as illustrated in Figure HI neither of these predictive relationships is as successful 
at predicting M, for elliptical galaxies as the M,{a) fit. Given that the bulge sample is 
statistically weaker and less stable (more subject to variation with fitting selections), as 
previously discussed, these two fits will not be further considered. The M,{Re, a) fit does not 
indicate M, oc a'^Re, the bulge mass surrogate, for any galaxy sample, nor does the M,{Re, a) 
fit appear superior to M,{a) for elliptical galaxies either from the fits, or from examining 
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the M,((j) residuals versus Re (Figure [3]). This is in contrast to the lMarconi fc HuntI (120031 ) 
argument, based on an examination of the residuals, that M,{Re,a) is superior to M,(cr) 
for all galaxies. Furthermore, neither M,{Re) nor M,{I(.) are better predictors of M,; for all 
galaxy samples they are worse (at a > 1-sigma level) than M,(cr). 

The three-parameter predictor, M,(/e, -Re, c"), is better (at the 1-sigma level) for bulges, 
and minimally better for ellipticals and the combined sample, when compared with M,(cr). 
The coefficients of the fit (Table [7]), however, vary strongly with galaxy sample; there is no 
clear relationship applicable for all galaxies. The bulge sample indicates a fit with minimal 
le and Re and a high a slope; this is a poor predictor for elliptical galaxies, as illustrated 
in Figure HI However, in both the elliptical and combined samples, the velocity dispersion 
is no longer the dominant term as it was in the M,(Je,cr) and M,{Re,a) projections of 
the multivariate surface. For elliptical galaxies where M, oc o'~^''^^I^'^^Rl'^^ , the preferred 
projection appears to be M,{Ie,Re) instead of M,((t). 

The M,{Ie,Re) fit results in a fit which is at least as (statistically) good as M,(cr). 
The form of the fit, for elliptical galaxies, (M, oc Jg'^^i?g'^^) is similar to the fit for the 
combined sample (M, oc /g'^^-R^'^^); both suggest M, oc leRl- The bulge fit predicts a 
different relationship, M, oc /g'^^i?g'^^, with a lower dependence on both variables. 

Examining the suggested M,{I^Rl) fit (Table [8|, Figure [5]), the scatter is minimally 
better than M,(cr) for all samples. The lower bulge-only slope is, within the stated un- 
certainties, consistent with the higher, eUiptical galaxy slope (0.76 versus 0.91-0.95); this 
behavior is reminiscent of the M,{a) results. This correlation is also physically significant; 
the SMBH mass is proportional to the energy (E) per T^, at the bulge effective radius: 
I^Rl oc (given E oc M^/Re and Le oc hRl). 

The correlation of M, with energy, M, oc T^I^R^ illuminates a problem with the T^pc 
mass-to-light ratios. The Tbpc values (see Section 12. 2p result in M,(T^/gi?g) fits which 
are worse (statistically) than M,(o"). This is illustrated in Table [8] as the Tf^^I^Rl pre- 
dictor. However, when th ese Thpc values a re rep laced with the (SAURON integral-field- 



spectroscopy) values from ICappellari et al.l (120061 ) . Tcap, when available, the fit improves 
significantly, as illustrated in Fi gure E] and in T a ble [5] (Tl^^^I^Rl). For the sample of 8 



elliptical galaxies for which the Cappellari et al. ( 20061 ) mass-to-light ratios are available 



(Cap8 \ M,{T'^I'jRl) is minirn ally better than both M,{I^Rl) and M,{a). It is anticipated 
that if Cappellari et al. ( 20061 ) (Tcap) values were available for the full sample of galaxies. 



the improvement in scatter would be significant; possibly enough to show a clear preference 



for M ,(E) over M,(cr). It is also interesting to note that although utilizing ICappellari et al. 



(120061 ) (Tcap) values improved the quality of the fit, the actual slope of the fit changed little 



or not at all. The elliptical and combined samples show a slope of 0.60-0.65, and the lower 
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bulge slope (0.47) is consistent with this range within the stated uncertainties. Figure [7] 
illustrates the higher, elliptical slope is a better descriptor of galaxies as a whole. Allowing 
for a multivariate fit M,{E, R^), tabulated in the bottom panel of Tabled does not improve 
the quality of the fits. 



Finally, a commonly proposed predictor of SMBH mass is luminosity. iMarconi fc Hunt 



mom determined that log (M./Mq) = (8.21 ± 0.07) + (1.13 ± 0. 12) {log {Lk, bulge) - 10.9) 
with an rms scatter of 0.31 (and similar slope and scatter in B, J and H). Using leRl as a 
surrogate for luminosity (Table [HI Figure [H]) it is clear that for ellipticals, M^^I^Rl) is inferior 
to M,{a) as a SMBH mass predictor. For bulges 'luminosity' is minimally better, although 
(see Figure [8]) there is still visible scatter about the fit if non-fitting sample galaxies are 
considered. For all samples the slope is roughly consistent with 1.0, as in the literature; the 
combined slope is in fact identical to the iMarconi fc HuntI (120031 ) K-band result with almost 
the same uncertainty. 



4.4. Gravitational Binding Energy 

The multivariate fitting [M,{Ie, Re,cr)] suggests that gravitational binding energy is a 
comparable, or better, predictor of SMBH mass than bulge velocity dispersion. The more 
formal fit of M,{Eg), where binding energy (Eg) is calculated utilizing previously discussed 
SB profiles (Section [2]) and geometrical assumptions (Section [3]), with no contribution from 
dark halos and disks, supports this conclusion. Throughout this discussion, and the accom- 
panying figures and tables, the term gravitational binding energy (Eg) is defined such that 
Eg = —PE, where PE is the gravitational potential energy. 

The most notable characteristic of the M,{Eg) fit is the stability of the slope with fitting 
algorithm, weightcent, and even with galaxy morphology. The best-fit slope is approximately 
0.6 (Table[9]) for the 15E, Cap8, and 23gal samples, both using the Isqfit and the robust fitting 
algorithms. This slope, 0.6, was also determined in the previous section using the crude 
approximation {T'^I^R^) for gravitational binding energy; the method for calculating the 
binding energy for the galaxy has no significant impact on the derived slope of the predictive 
relationship. The variation in the slope based on the weightcent selection, as illustrated by 
the ttmax and amin values in Table [H is minimal. Even the addition of uncertainty in the Eg 
value, = 0.3 (100%) (Table [TU]) using the fitexy algorithm, results in a slope which is, at 
most, 0.03-0.04 higher and within the stated uncertainties for the Isqfit (~ 0.6) slope. 

The bootstrap algorithm results, implemented as for the M,{a) predictor, further illus- 
trate the stability of the slope. For the 15E sample, the mean and the standard deviation 
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of the slope are identical to the original (single-run Isqfit) results when using the bootstrap 
Isqfit algorithm (Table [TT]) . The remaining samples {8B, Cap8, 23gal) show only slightly 
lower slopes and slightly higher uncertainties in the Isqfit bootstrap slopes, when compared 
with the original slopes. The large differences between the original and bootstrap slopes 
(and uncertainties) for the 8-galaxy samples {8B, CapS) which were present in the M,{a) 
predictor are not present here; the Eg predictor is more stable even for small sample sizes. 
The robust bootstrap (Table [T2!) is also generally consistent with the original results, for the 
elliptical and combined samples, further strengthening the argument for M, cx E^-^ as a very 
stable and reliable predictor. 

The spiral/lenticular bulge-only fitting sample {8B) slope is consistently lower than for 
the elliptical samples for the least-squares fitting; this was also seen for the M,{a) predictor. 
This lower slope is given little credence, however, based on the previously mentioned weak- 
nesses in the bulge sample, the fact that the Isqfit slope of 0.47 is consistent with 0.6 within 
the stated uncertainties, and the fact that the slope is not stable with fitting selection. The 
robust fitting algorithm predicts a higher slope (0.66), while the robust bootstrap produces a 
slope of 0.43, consistent with the Isqfit bootstrap results. This discrepancy is indicative that 
the bulge sample is a much weaker (less consistent) predictor in comparison with elliptical 
and combined galaxy samples. The preferability of the ~ 0.6 fitting slope for galaxies, in 
general, is further illustrated in the top panel of Figure [2] and in the correlation between the 
5i?-fit-residuals and Eg for the elliptical galaxy sample (Figure [9]). 

The M,{Eg) predictor is not only at least as stable a predictor as M,(cr) with mor- 
phology, fitting algorithm, and weightcent, but it is also comparable in a statistical sense. 
The scatter is statistically equivalent to M,{a) for all samples, and minimally better for the 
elliptical {15E, Cap8) samples. The comparable nature of these two predictors, particularly 
for ellipticals, is visually illustrated in Figure [2], in which M,{Eg) is given along the top and 
M,(ct) along the bottom. This similarity in predictor-quality is remarkable given that a is 
based on direct observations, while Eg is dependent upon a series of geometrical assumptions 
and the combination of (often discrepant) data; an error in any of these assumptions or data 
combinations will increase the scatter in a manner not present for the simple a measurement. 
As illustrated in the previous section, improved T values for all 23gal sample galaxies will 
likely reduce the scatter. There is no indication for the necessity of either a second parameter 
(Re) or for a quadratic term in the residuals (Figure [9]), as further discussed in Appendix iBl 



Th e replacement of the SMBH mass for NGC 224 with the higher value of 1.4xlO'^M(r 



given m 







Bender et al.l (120051 ). and the replacement of the 140 kms~^ velocity dispersion, for 
NGC 2787, with the higher value of 200 kms~^, does not change the conclusions drawn from 
the fits, as shown in Table [91 The slope of both the 23gal least-squared {Isqfit) fit and the 



- 23 - 



23gal robust fit are unchanged. The slope of both the 8B Isqfit and robust fits are only 
marginally higher; the least-squared fit is raised to 0.49 ± 0.14 from 0.47 ± 0.15, and the 8B 
robust fit slope is raised to 0.71 from 0.66. The scatter and requisite intrinsic uncertainty 
in logM, {eyin) is smaller for all but the 8B robust fit (for which the scatter is only slightly 
higher), which, when combined with the increase in scatter for the M,{a) fits using these 
revised values (see Table [5]) reinforces the relative success of the gravitational binding energy 
in accurately predicting SMBH masses. 



4.5. Gravitational Potential 

The predictive power of the gravitational potential, $(-R), is also explored. This quantity 
is determined as a by-product of the binding energy calculation (see equation [3]). The 
gravitational potential is evaluated as a SMBH mass predictor, M,($(i?)), at four radii: 20 
pc, 100 pc. Re, and -Re/8 (Table [T^. 

Examination of the fits at fixed radii (i.e. 20 pc and 100 pc) indicates that there is 
no single predictive relationship for the black hole mass which is applicable to all radii 
or to all morphological classes. The strength of the fit decreases as radius increases for 
the elliptical samples {15E, CapS), but not for the combined or spiral/lenticular galaxy 
samples. For all samples (but most notably for the elliptical samples) the slope decreases 
with increasing radius. The variations in slope with weightcent, as given by amax and amin 
(Table [T3!) . are within the stated uncertainties for each sample, but the range in these values 
indicates that the predictor is not as stable as the gravitational binding energy or bulge 
velocity dispersion. Furthermore, the elliptical sample slopes are irreconcilably higher than 
the bulge sample slope for both the 100 pc and the 20 pc radii fits; these are two separate 
populations. The combined sample slope is intermediate between the elliptical-only and 
the spiral/lenticular bulge-only slopes, but it is a poor predictor; it is worse than velocity 
dispersion at a statistically significant level. 

The M,($(ni?e)) predictor, evaluated at Re/8 and Re, does not alleviate the problem 
of an irreconcilable slope between the spiral/lenticular bulge and the elliptical samples. The 
combined sample slope is a compromise between the high (~2) slope from the elliptical 
samples and the low (~1) slope from the bulge sample. The result of fitting these two 
disparate populations simultaneously is a predictive fit which is (statistically) worse than 
velocity dispersion. 

The M,($(?7.i?e)) predictor does, however, exhibit several features common to all mor- 
phological samples. First, as the radius decreases, the slope of the predictor increases, but 
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only slightly, and to a value within the original stated uncertainties. Second, as the radius 
decreases, the scatter decreases, meaning that the strength of the predictor increases. Third, 
the robust fitting slopes are higher than the least-squares fitting slopes, irreconcilably so for 
the spiral/lenticular {8B) samples. Thus, this predictor is less stable than the energy or 
velocity dispersion predictors. 

The M,($(i?e/8)) predictor for the elliptical samples is minimally better than the (sta- 
tistically equivalent) M,{a) predictor; however, it is a much less stable predictor. The 
statistical equivalence is expected since $^ ~ a^; in essence these two predictors are fun- 
damentally related. The lack of stability is illustrated by several characteristics. First, the 
robust slopes for the 15E and Cap8 samples do not agree with each other. Second, there is 
a noticeable range in slope encompassed by amm and amax for the Isqfit results. Third, when 
uncertainty in $(i?e/8) is included (Table [TUl). there is often no viable eyin selection to pro- 
duce Xr = 1-0, and higher slopes result when a fit is possible. Finally, the bootstrap fitting 
(Tables [11] - [12]) shows discrepancies between the mean and the median slope, particularly 
for the Cap8 sample. 

The M,($(i?e/8)) predictor for the spiral/lenticular bulge sample is not an improvement 
on M,(cr). First, based on the scatter, it is minimally worse for the Isqfit predictor, and 
statistically significantly worse for the robust predictor. Second, like the elliptical sample 
there is no slope stability between the robust predictor slope and the Isqfit slope. Third, 
there is a discrepancy between the robust bootstrap mean and median slope. Finally, when 
large uncertainty in is added, the slope of the predictor becomes irreconcilably higher. 

The large variations in both the original fits and in the bootstrap results with fitting 
algorithm and weightcent selection, combined with the high slope for ellipticals (> 2.0) and 
the low slope for the bulges (0.98-1.37), as illustrated in Figure [TUl are not reconcilable 
within the estimated uncertainties, and argue against gravitational potential being a strong 
or stable predictor of black hole mass for most galaxies. It is, at best, an adequate predictor 
to use for elliptical galaxies only. 

The fit residuals plotted against both $(i?e/8) (Figure [TT] left) and Re/S (Figure [TT] 
right) show neither a clear need for a quadratic, nor the need for an additional parameter, 
as confirmed in Appendix [B] However, the correlation of the iJ-E-fit-residuals with $(/2e/8) 
for the spiral/lenticular galaxies and the 55-fit-residuals with $(i?e/8) for the elliptical 
galaxies reconfirms the assertion that these are different populations with irreconcilable 
slopes. The slight ^%a/-fit-residual correlation with for the elliptical galaxies is a result of 
the 'compromise' combined {23gal) slope; the elliptical-only-fit-residuals (top panel) do not 
show this correlation. It is also notable (Appendix [B]) that the M,[$(i?e/8), -Re] multivariate 
fit for spiral/lenticular galaxies suggests a property akin to mass may be slightly preferred 
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to potential alone. 



4.6. Bulge Mass 

The galaxy bulge mass predictor has been found by previous authors to correlate 
(roughly) linearly with SMBH mass. In previous studies it was found that log (M./M©) = 
(8.28 ± 0.06 ) + (0.96 ± 0.07) (log Mw „. - 10.9) where M^^ige = 3Rea^^G-\ with an rms scat- 



ter of 0.25, ( Marconi &: Hunj 2003), and, for a 90 ga laxy sample of both active and inactive 



galaxies, M. oc M^^^f^-^^ dMcLure fc DunlopI 120021 ). Using more detailed r nass- modeling 



(eithe r the Jeans equation or alternative dynamical models in the literature), iHaring fc Rix 
J2OO4I ) found log {M,/Mq) = (8.2O±O.lO) + (1.12±O.O6)log(M6„/ge/lO^^M0), with a scatter 



of 0.3 dex, using a bisector linear regression fit (and an error in Mbuige of 0.18 dex). When 
fitexy, the algorithm used here, is applied, they quote a slope of 1.32 ± 0.17. 

The galaxy mass explored here as a predictor of SMBH mass is the bulge stellar mass 
enclosed within a sphere of radius R, M{R), derived from the SB profile (Section 12. ip and 
the previously mentioned geometrical and physical (i.e. no dark halo or disk contribution) 
assumptions (Section [3]). As with the gravitational potential, this enclosed mass is examined 
both at fixed radii (1 kpc and 10 kpc) and as a function of the galaxy size [Re and lORe). 
Of these four radii, lORe is the best substitute for the bulge mass; there should be very little 
galaxy mass beyond a few Re, and so all of the mass from the bulge should be enclosed by 
the extreme limit of 10-Re- 

The predictors based on the mass enclosed by the fixed radii (i.e. 1 kpc and 10 kpc) 
are not an improvement on the bulge velocity dispersion predictor. For the combined galaxy 
{23gal) sample, the M{lkpc) predictor is statistically worse than M,{a). The M{lkpc) fit 
for the remaining samples, and the M{10kpc) fit for all samples, exhibit scatter which is 
statistically equivalent to, but minimally worse than, the M,(cr) predictor (Table [T^ . The 
comparison of the M{10kpc) and M{lkpc) predictors indicates that the M{10kpc) predictor 
is better (it has minimally lower scatter for all but the 15E sample and a more stable slope 
with morphological selection) than the M{lkpc) predictor. The M{10kpc) predictor slopes 
are stable (at 0.9) with both morphological sample and with weightcent selection. Given the 
consistency of this slope and the fact that statistically it is equivalent to the M,{a) predictor, 
it is adequate as a SMBH mass predictor, but it is not a significant improvement over either 
velocity dispersion or gravitational binding energy. 

The M{Re) and M{10Re) predictive fits probe a more physically uniform region of the 
galaxies than the fixed-radius fits, but do not produce a predictor which is an improvement 
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on velocity dispersion. The scatter is never lower than for the velocity dispersion predictor. 
The M{Rf.) 15E and 23gal sample predictors have statistically significantly worse fits, while 
the remaining fits are generally statistically equivalent to M,(cr). The M(lO-Re) predictor 
fits always have lower scatter than for the M{Rf,) fit. Given that the M{Re) predictor is 
weaker than the M{10Re) predictor, it will not be discussed further, other than to note that 
the slopes are consistent with 0.8 for all morphologies and weightcent selections. 

The slope of the M(10i?e) predictor is very consistent (~ 0.8), for the 15E, Cap8 and 
23gal fitting samples, with morphology, fitting algorithm and weightcent selection. The 
robust slopes are slightly lower, 0.77 in comparison with 0.81-0.83, but well-within the stated 
uncertainties. The bootstrap fitting (Tables [11] - [12]) further illustrates the stability of the 
0.8 slope for these three sam ples. In general, the fitted rela tionships are roughly lin ear and 
the slopes are consistent with lMcLure k. Dunlop (I2OO2I) and iMarconi &: HuntI (j2003l ). within 



Haring k Rixl (120041 ) . likely due to slight 



the stated uncertainties, and slightly lower than 
differences in the mass determination, sample selection and fitting. 

Although it exhibits a consistent slope, the M,(M(10-Re) predictor is not as stable overall 
when compared with predictors such as M,{Eg). First, the robust fits are statistically worse 
than the velocity dispersion predictor for both the elliptical {15E) and the combined {23gal) 
galaxy samples. Second, the addition of = 0.3 (100%) (Table [10]) results in a noticeable 
slope increase of 0.12 for the 23gal sample and an increase of 0.26 for the 8B sample, both 
outside of the original stated uncertainties. Even for the full-elliptical sample {15E), the slope 
is increased by 0.10, at the limits of the stated uncertainties. While these may be extreme 
limits, it exhibits a lack of the stability in slope, in comparison with the gravitational binding 
energy prediction. Third, the spiral/lenticular bulge sample {8B) has a slightly higher Isqfit 
slope, 0.88, well within the stated uncertainties, but shows inconsistencies by exhibiting a 
noticeably higher slope with the robust fitting algorithm (1.23), and a noticeably lower slope 
(0.75) with the Isqfit bootstrap fitting. It does center on the more common 0.8 slope for the 
robust bootstrap fitting. 

The fits for M, versus (M(10i?e)), shown in Figure [12], illustrate the similarity in the 
slope for all of our galaxy samples. The residuals of the fits (Figure [12]) , plotted against 
(M(lO-Re)) and 10i?e show neither a strong indication for a quadratic, nor for the addition 
of another parameter for the bulge and combined samples; this is reconfirmed in Appendix [Bl 
The addition of Re does improve the scatter in the elliptical samples: the best fit is M, oc 
(M(10-Re))^'^-R~^'^, an indication that pure elliptical galaxy bulge mass is not the strongest 
predictor; rather a physical quantity closer to potential or binding energy is suggested. 
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5. DISCUSSION & SUMMARY 

The goal of these calculations was to determine (1) whether the addition of a second (or 
third) host-galaxy-bulge parameter to the bulge velocity dispersion predictor of SMBH mass 
would result in an improved predictor (smaller residuals) and (2) if an alternative predictor 
(such as gravitational binding energy or gravitational potential) is equivalent to, or superior 
to, the well-accepted velocity dispersion predicto r. These calcula t ions w ere undertaken for 



four (sub)samples of galaxies, all taken from the iTremaine et al.l (120021 ) 31-galaxy sample: 
a sample of 15 elliptical galaxies {15E), a sample of 8 spiral/lenticular bulges {8B), the 
23 galaxy combination of these two s amples (23 gal), and fin ally, a subsample of 8 elliptical 



galaxies ( Cap8 ) with T available f romlCappellari et al 



predictor, using dispersions from ITremaine et al. j2002f L was evaluated for each of these 4 



2006h . The bulge velocity dispersion 



samples to provide a reference point against which the other predictors were compared. 
The remaining predictors were based on SB profiles and galaxy parameters taken from the 
literature, along with the assumption of an oblate spheroidal geometry with no disk or dark 
halo contributions. Thus, the comparisons between the predictors are based on identical 
samples and the same data set and assumptions, allowing for a direct comparison of 
the predictor strengths and weaknesses. 

In answer to the first question posed, no additional parameter, when combined with a, 
was found to produce a better predictor than M,(cr) for the elliptical and combined galaxy 
samples; however, the only galaxy parameters examined in a multivariate fit were le and Re. 
Multivariate fits of M,{Ie,(j) and M,{Re,a) were not better than M,{a) for the 15E and 
23gal samples, and had negligible coefficients for I^. and Re- The combination of all three 
parameters (M,(/e, ReyCr)) for these two samples yielded the suggestion that an alternative fit 
depending only on M,{Ie, Re) was warranted. The M,{Ie, Re) predictor shows an indication, 
for ellipticals and the combined sample, that M, oc [Jfi?g]^ oc E^, where 'E' is the binding 
energy of the galaxy bulge. 

The M,((t) fit was also examined as function of morphology. While the predictive-fit 
slopes were different (the spiral/lenticular sample fit exhibits a lower slope than the elliptical 
samples), they were statistically reconcilable with each other. However, the M,{a) residuals 
for the spiral/lenticular galaxies do correlate with Rg, suggesting the nec essity of a multi- 



variat e fit, a phenomenon which is not suggested for the elliptical galaxies. iMarconi fc Hunt 



(120031 ) likewise identified a correlation between the M,{a) residuals and i?e, but did not dif- 
ferentiate based on galaxy morphology, and attributed the correlation to all galaxies, rather 
than to just the spiral/lenticular galaxies which visibly show the correlation. 

In answer to the second question posed, an alternative predictor to M,((t), gravitational 
binding energy, was found to be as (statistically) good as the a-predictor for all samples. 
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Considering that it is equivalent, even with the simphfying assumptions and data constraints, 
this r aay be an indica t ion t hat gravitational binding energy is a fundamental predictor. 



Using ICappellari et al.l (120061 ) mass-to-light ratios for available galaxies, the scatter becomes 
noticeably lower. This implies that if such mass-to-light ratios, or other improvements in 
the data, were available for all galaxies, the scatter could be further reduced, potentially 
making this predictor (statistically) superior to bulge velocity dispersion. It may further be 
hypothesized that since the modeled-Eg predictor, with its many simplifying assumptions and 
data imperfections is as good as M,(cr), that if the real, physical- E'^ were known precisely for 
each galaxy, as the velocity dispersion is known, the Eg predictor would become (statistically) 
superior to velocity dispersion. Given that replacing the surrogate-energy (T^/gi?g) with the 
more formally calculated Eg decreased the scatter but did not have an impact on the slope 
of the Eg predictor, it is anticipated that replacing the modeled-i?^ with the real-Eg would 
not alter the fundamental relationship (and slope), but only the predictive strength. 

The slope of the predictor, M, oc -E^'^, is remarkably stable. It shows minimal variation 
with changes in fitting algorithm (least-squares versus robust), least-square centering of M,, 
least-square weighting selection (both M, and Eg), or method of calculation of Eg (crude 
[T^/gi?g] versus formally calculated). Calculating the energy as M(10i?e)o'^ or TlgRj^a'^ 
produces the same slope, 0.6, as well. Even utilizing the T;,pc-values does not affect the 
slope. The bootstrap (with both least-squares and robust algorithms) further reinforces 
the slope stability. There is, however, some variation with morphology. The bulges {8B) 
have a lower (0.47) least-squares slope [similar to M(cr)], but this is equally likely to reflect a 
statistical weakness in the sample as a true morphologically-based difference in the predictive 
relationship. 

In general, the spiral/lenticular bulge {8B) sample predictors are inferior to the 15E or 
23gal predictors. First, the sample size is (statistically) significantly smaller; the 8B sam- 
ple is ~ one-half the size of the 15E sample and ~ one-third the size of the 23gal sample. 
Second, for the 8B sample, in contrast with the 15E sample, observational-uncertainty-only 
weighed least-squares fits give very different results from the unweighted fits, suggesting 
errors in M, (or the associated observed uncertainties). Eliminating additional 8B sample 
galaxies to alleviate the problem is not feasible as this would further weaken the sample in a 
statistical sense. Third, other than the velocity dispersion predictor, all of the predictors are 
dependent on a (literature-based) bulge-disk decomposition with the assumption that the 
galaxy has only two components. These added structural and geometrical assumptions for 
the spiral/lenticular galaxies may have a weakening effect on the predictors; false assump- 
tions or over-simplifications could make a predictor appear weak for the spiral/lenticular 
bulge sample, when in truth, it is the structural simplifications which are at fault, not the 
fundamental predictor. In this sense, the assumption-free, directly observable velocity dis- 
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persion may be a more "reliable" predictor of SMBH mass than any o f the other predictors. 
This is similar to the conclusion reached by lMcLure &: DunlopI (120021 ) . 



The calculations of the gravitational binding energy and the gravitational potential 
omitted the contribution from the spiral/lenticular disk. Examining Figure [2] (for the grav- 
itational binding energy) and Figure [TO] (for the gravitational potential), it is obvious that 
while the spiral/lenticular bulge {8B) and elliptical {15E) populations do not have the same 
slope (the bulges prefer a slightly shallower slope), they are occupying the same region of 
parameter-space; the bulge-population is not offset from the elliptical-population, and the 
two predictive relations cross in the center of the fitting range. The inclusion of the spi- 
ral/lenticular disk would increase the galaxy mass and energy and likely result in an offset 
between the elliptical and spiral/lenticular galaxy populations; however, the formal calcula- 
tion of the impact of including the disk-contribution is beyond this scope of this paper. 

For the observational h ost galaxy-SMBH cor relations in the literature, includin g the 



log- linear ve locity dispersion (ITremaine et al. 



20021). luminos i tv flMarconi fc Huntll2003r) . and 



bulge mass jMcLure fc DunlopI I2OO2I : iMarconi fc Huntlbood : iHaring &: Rixll2004h predictive 



relationships and the log-quadratic velocity dispersion predictor (lWyithell2006l ). the fits ob 



tained here are generally consistent with the published results. The log-linear velocity dis- 



persion predi ctor determinec 



than that in Tremaine et al 



lere for the full galaxy sample exhibits a slightly lower slope 
(120021). bu t is fu lly consistent within the stated uncertainties 



and with statements in iTremaine et al.l (|2002| ) regarding the slope as a function of sample 
selection. The luminosity was not investigated here, but the luminosity-surrogate (/e-Rg) P^^" 
dictor was found to have a r oughly linear rela t ionsh ip with SMBH mass, and to exhibit the 
same slope as determined in iMarconi fc HuntI (120031 ) for the luminosity. The bulge mass, as 
estimated here by summing all of the visible stellar bulge light under the assumptions of an 
oblate spheroidal geometry for the bulge and a constant mass-to-light ra tio, is found to cor- 
relate roug hly linearly with SMBH mass, as also seen in the literature by iMcLure fc Dunlop 
(I2OO2I ) and IMarconi &: HuntI (120031 ). The bulge masses, de termined using the Je ans equation 
and other dynamical models from the literature, used by iHaring &: Rixl (120041 ) led them to 
find a slightly higher slope than we determine here; however, we do not believe the differ- 
ence to be sign ificant. Finally , we find log-quadratic fits which are consistent with the fits 
determined by IWyithd (120061 ). However, using our indicato r of stat i stical superiority (the 
Snedecor F-test instead of the Bayseian analysis employed in IWvithel J2006h ) we do not find 
these fits to be better than the log-linear fits. 



Our findings are also generally consistent with the theoretical results from lHopkins et al. 



(l2007l ). who probed predictors of SMBH mass, and the SMBH mass fundamental plane, using 
major galaxy merger simulations. They find, in their simulations, that M, oc [M^.a^]'^'^^^'^'^^, 
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with a smaller (simulated) scatter than for their M,(cr) or M,(M*) relations, where 
is the stellar mass. This result is consistent, within the stated uncertainties, with our 
observationally-based determination that M, oc ii;0-62±o.o6 f^j, elliptical galaxy popu- 
lation [15E), and confirms a 'tilte d' relationship betwee n supermassive black hole mass and 
gravitational binding energy. The iHopkins et al.l (120071 ) simulations also find the existence 
of a supermassive black hole fundamental plane, M, oc o-2.9o±o.38^o.54±o.ii^ ^^^^ claim this 
to be better than any single-variate predictor of SMBH mass, as determined examining 
a, M^, Mdyn, or R^, at a greater than 3-sigma level. Their stated scatter for this rela- 
tionship (0.21), however, is identical to that given for their gravitational binding energy 
predictor; we estimate the best predictor using a comparison of the scatter, and so under 
our statistical methodology their gravitational binding energy predictor would be considered 
equal to this fundamental plane predictor. In examining our (weaker) spiral/lenticular {8B) 
sample, we find M, oc cr2.90±o.52^o.94±o.33 to be better than M ,(a) at a statistically signif- 
icant level, which is consistent with the iHopkins et al.l (120071 ) fundamental plane, within 
the stated uncertainties. However, for our sample containing only ellipticals {15E), we find 
M, oc (j3 '^3ifco.7i^o.o5iii0.24 .^j^jg r esult has a substantially different dependence on than that 



found by IHopkins et al.l (120071 ) , and we do not identify this as being an improvement on the 
M,(a) predictor. Th e reasons for this discrepancy, in ascertaining the best predictor, with 
Hopkins et al.l (120071 ) likely stem, in part, from the difference between the modeled/simulated 
galaxy parameters and the observationally-determined galaxy parameters, which inherently 
possess additional scatter that may obscure a fundamentally multivariate relationship. 

The main conclusions of this study can be summarized as follows. 



Gravitational binding energy (M, oc E[ 
dispersion. 



0.6^1 

9 ' 



is as strong a predictor of M, as velocity 



(a) The scatter is statistically equivalent for both M,{Eg) and M,(cr). 

(b) The slope is very stable and does not vary with data-weighting selections, fitting 
algorithm applied (least-squares or robust), or with which sample galaxies are 
selected for inclusion in the fitting (bootstrap). The slope resulting from the full 
data set (spiral/lenticular bulges and ellipticals) is the same as that for elliptical 
galaxies only. 

(c) Energy is suggested by multivariate M,(Je, Re, <j) fitting for ellipticals. 

2. The spiral/lenticular bulges and elliptical galaxies may be different populations. 

(a) M,((t) residuals correlate with Re for spiral/lenticular g alaxie s, but not for el- 
lipticals. This is in slight contrast to iMarconi &: HuntI (120031 ) who identified a 
correlation but did not differentiate based on morphology. 
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(b) Multivariate fitting suggests different physical predictors for bulges (a mass-like 
quantity) and ellipticals (an energy-like quantity). 

(c) Spiral/lenticular bulges have lower (but statistically reconcilable) slopes for M,{Eg) 
and M,((j), and substantially different and irreconcilable slopes for the gravita- 
tional potential predictors. Some of these differences may stem from a statistically 
weaker spiral/lenticular bulge population sample when compared with the ellip- 
tical sample. 

3. Bulge mass is an adequate predi c tor o f M,, with a rough l y line ar slo pe, confirming 
the re sults of iMcLure fc DunlopI j^OoJ, iMarconi fc HuntI (EooJ, and iHaring fc Rix 
( 120041 ). However, the bulge mass predictor is not as stable with variations in fit- 
ting algorithm (least squares versus robust) or with data-weighting selections as either 
M,{Eg) or M,((t). Furthermore, the multivariate fit for elliptical galaxies suggests a 
combination of mass and Re, closer to energy, is preferred over mass. 

4. Gravitational potential is an adequate predictor of M, for elliptical galaxies. However, 
it provides no new insight over the comparable M,(ct) predictor and it is an inferior 
predictor for spiral/lenticular bulges when compared with M,{Eg) and M,(cr). Fur- 
thermore, the spiral/lenticular bulge and elliptical populations exhibit irreconcilably 
different slopes, and there is a slight indication from multivariate fitting that a mass- 
like quantity is preferred for the spiral/lenticular bulges. The slope is not very stable 
with variations in fitting algorithm (least squares versus robust) or data-weighting 
selections for any of the fitting samples. 

5. 'Luminosity' (/e-Rg) is an adequate predictor of M, with a roughly linear slope, as 
seen in the literature (IMarconi &: HuntI l2003l ). It is statistically equivalent to M,(o") 
for the spiral/lenticular bulge sample, and for the combined galaxy sample, but it is 
statistically worse for the elliptical galaxy sample. 

6. Improved values of the mass-to-light ratio would improve the scatter in the M,{Eg) fit. 



7. There is no strong evide nce (see Appen dix IB.1|) for the need for a log-quadratic fit- 
ting function, preferred by lWyithd (120061 ) . for any of the parameters examined (velocity 
dispersion, gravitational binding energy, gravitational potential, or bulge mass). Us- 
ing the criteria specified in this paper to select the best form of the predictors (the 
Snedecor F-test), the log-quadratic predictors are statistically equivalent to the log- 
linear predictors. Therefore, the simpler, log-linear predictors are pref erred. However , 
the form of the log-quadratic predictors are compatible with those in IWyithd (120061 ) . 
and it is possible that with a sample that includes more high and low mass galaxies. 
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that the necessity for the log-quadratic predictor would become more apparent using 
our selection criteria. 
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A. SENSITIVITY TO VARIATIONS IN AXIS RATIO, INCLINATION, 

AND MASS-TO-LIGHT RATIO 

In modeling the galaxy as an oblate spheroid in order to calculate the gravitational 
binding energy, the gravitational potential, and the bulge mass, it is crucial to have accurate 
estimates of the input parameters for each galaxy (see Section [2^ . In the literature, however, 
there is often a wide range of published values for a given galaxy, particularly for the axis 
ratio, inclination and mass-to-light ratio. In this section, the effects of having selected the 
"wrong" value for each of these input parameters, is explored. In addition, those galaxies 
with the largest likelihood of having such an incorrect value selection are identified. 

The intrinsic axis ratio (p) is used to project the observed major axis surface brightness 
(mass) profile to a three-dimensional description of the galaxy (see Sections [2] - [3]) . The 
intrinsic axis ratio depends on both the observed inclination {6) of the galaxy and on the 
observed axis ratio (q). The effect on the calculated gravitational binding energy of varying 
each of these observed quantities over the full range of physically allowable values, 0.0 < 
q,p < 1.0 and 0.0° < 9 < 90.0°, is illustrated in Figure [T^ for three elliptical galaxies: NGC 
221 (the smallest galaxy), NGC 4486 (the largest galaxy and one of the most inclined) and 
NGC 6251 (an intermediate galaxy). In this plot q/qtase (bottom panel) and 9/9base (top 
panel), where the 'base' value is the value determined in Sections [2] - [3l is plotted against 
Eg/Eg^^^^. Only one of the observed quantities is varied at a time while the other is left fixed 
to the value given in Section [21 the two observed quantities are never simultaneously varied. 

The inclination (9) is difficult to determine for the diskless elliptical galaxies, which are 
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usually taken to be at, or near, edge on, while for the spiral/lenticular galaxies the disks 
provide a better constraint, but can still be affected by assumptions about the disk thickness 
or the presence of disk warping or non-round disks. Those galaxies with a literature-value 
10% larger or smaller than the best value used in Section [2] are given in Table [T3 (It is 
possible, throughout this section, that there exists an unidentified literature source which 
would cite a still more discrepant value, or a discrepant value for one of the other sample 
galaxies.) Based on Table [151 Eg/Eg^^^^ = 1.1 for NGC 221 {9 = 50°), Eg/Eg^^^^ = 0.98 for 
NGC 4486 {9 = 51°; for 9 = 90° it would be 0.95), and for NGC 6251 the variation would 
be negligible. Inspecting Figure [HI it is evident that incorrect inclinations, even those off 
by a substantial amount, will have a minimal effect on the calculations. 

The observed axis ratio (q) also exhibits variation in the literature-quoted values which 
may reflect an axis ratio (ellipticity) that is dependent on the isophotal fitting algorithm, a 
strong color-dependence of the shape (and ellipticity) of the galaxy, or a physical structure 
in the galaxy which renders the axis ratio not constant with radius (in violation of the 
simplifying modeling assumptions adopted here). The galaxies for which there is a bulge 
q- value in the literature differing from the best value by 10% or more are given in Table [T6l 
The plots (Figure [T3J bottom panel) all exhibit a quadratic relationship between the ratio 
of the axis ratios {q/qtase) and the ratio of energies {Eg/ Eg^^^J. Based on Table [16], the 
range of observed Eg/Eg^^^^ values would be 0.8-1.3 for NGC 221, 0.9-1.1 for NGC 4486 (the 
stated lower limit of q=0.65 is unphysical at an inclination of 42°), and 0.97-1.2 for NGC 
6251. Looking at the full range of Eg/Eg^^^^ for these representative galaxies, the binding 
energy could physically vary by an order of magnitude, but only in the case of an axis ratio 
which is in error by a factor of two, which is not seen for these galaxies and which is unlikely 
for any of the galaxies based on Table [T6l For the majority of the galaxies, the variations in 
the energy, based on q, are at most 20-30%, which is accounted for in the modeling. 

Finally, the calculated gravitational binding energy is strongly dependent on the mass- 
to-light ratio adopted; E oc oc T^, so for an error in T such that T = xTtme, E = x^Etme- 
The values quoted in the literature show variations resulting from different estimates of the 
bulge luminosities, different methods for dynamically estimating the bulge mass, and possibly 
from a radial variation in the mass-to-light ratio. The least-well-determined T values, in the 
literature, tend to be for the spiral/lenticular bulges. It is possible to obtain an estimate of 
Mhuige from the rotation curve for these galaxies to produce an alternative T estimate, as 
was done for three of the galaxies (NGC 2787, 4459 and 4596), but for the majority of the 
galaxies, the bulge-only light profile, which is believed to be a source of scatter in and of 
itself, will result in an T which is, at best, a minimal improvement over the virial/dynamical 
mass determinations from the literature. In Table [TTl those galaxies for which there is a 
mass-to-light ratio in the literature > 10% different from the Tbpc value for the mass-to- 



-34- 



light ratio given in Table [2] are listed, and the most extreme mass-to-light ratios located in 
the literature are enumerated in order to probe the full range of possible variations which 
might exist in the computed energy. The ICappellari et al.l (120061 ) mass-to-light ratios are 
excluded from consideration in the construction of this table, although they are taken to be 
superior and are used for all galaxies where available (see Section I^l3|) . Table [T7I shows that 
for a few of the galaxies, if the extreme limits to the mass-to-light ratio are adopted, there 
could be a substantial variation in the computed gravitational binding energy. However, in 
including large uncertainties in the energy during the course of the fitting, it was found that 
this would not have a substantial effect on the derived M,{Eg) predictor (see Section 



B. LOG-QUADRATIC AND MULTIVARIABLE FITTING RESULTS 

B.l. Log-Quadratic Fits 



(jWyithell2006l ) suggests that both M,{a) and M,{Mbuige) are not log- linear functions, but 
are, in fact, log-quadratic expressions such that logM, = (8.05±0.06) + (4.2±0.37) log((T„) + 
(1.6 ± 1.3) [log (cT„)]2 (intrinsic scatter 0.275 ± 0.05) where cr„ = a/200kms-^ and logM, = 
(8.05 ± 0.1) + (1.15 ± 0.18) log (Mbn) + (0.12 ± 0.14) [log (Mft^)]^ (intrinsic scatter 0.41 ± 0.07) 
where Mf,n = Mbuige/^O^^MQ. This argument is based in part upon a rigorous statistical 
analysis of the data, and in part on the positive log- linear-fit residuals for both the smallest 
and the largest galaxies in the samples. 

For the sample of galaxies used here, the M,{a) residuals for fitting-sample galaxies, 
when plotted against a (Figure [3]), show no indication for a pattern of positive residuals 
with high- and low-cr galaxies in either the elliptical galaxy or in the spiral/lenticular galaxy 
subsamples. Likewise the fitting-sample galaxy M,{Eg) residuals vs Eg (Figure [9]), the 
M,{M{10Re)) residuals vs M(10i?e) (Figure [13]) and the ^%a/-fit M,{<!>{RJ8)) residuals 
vs $(i?e/8) (Figure [TTi) show no evidence for a pattern of larger positive residuals for the 
high and low mass galaxies. (The (?5-sample-fit A'/,($(i?e/8)) residuals vs $(i?e/8) for 
the elliptical galaxies and the iJE'-sample-fit M,($(i?e/8)) residuals vs $(-Re/8) for the 
spiral/lenticular galaxies exhibit inadequate fits; this is evidenced by linear correlations, but 
there is no indication of a quadratic). However, for all of these predictors, the three largest 
elliptical galaxies exhibit positive residuals, when compared with other sample-galaxies. 

When log-quadratic fits are performed for these four predictors (X= a, Eg, $(i?e/8), and 
M(10i?e)), using the multivariate fitting algorithm Isqfit, with normalized fitting parameters 
such that Xn = X/Xq (see Table [3]), there is no strong evidence for the necessity of a 
log-quadratic expression. The log-quadratic fits are tabulated in Table [181 where for the 
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function 

log M. = ai [log + as log X„ + d (Bl) 

the best fit coefficients (oj, d) are given for each of the morphologically-separated fitting 
samples {15E, Cap8, 8B, and 23gal), along with the requisite intrinsic uncertainty, the scatter, 
the maximum and minimum ai coefficient obtained by varying the weightcent selections, and 
the F-test ratio and its statistical significance (see Section H]), both when compared with the 
log-linear M,(cr) predictor and when compared with the log-linear M,{X) predictor. For 
every predictor and every fitting sample the log-linear and log-quadratic fits are statistically 
equivalent, using our statistical indicator, the Snedecor F-test. 

For the M,{a) predictor, the coefficient of the quadratic term (ai) varies strongly with 
morphological fitting sample and exhibits uncertainties larger than the coefficient itself. For 
the combined galaxy sample the best fit is logM, = (8.13 ± 0.08) -|- (3.90 -|- 0.37) log(T„ + 
(0.97 ±1.22) [l og (7nP (intrinsic uncertainty 0.25). These coefficients are consistent with the 



Wvithd (120061 ) result 



For the M,{Eg) predictor, the coefficient of the quadratic term for all of the fitting 
samples is minimal. For the combined sample the best fit is log M, = (8.06 ± 0.09) + (0.61 ± 
0.06) log E g ^+ {0.05 ±0.04) [log Eg ^]'^ (intrinsic uncertainty 0.27). The dominant term remains 
M, oc Eg'^, as previously seen in the log-linear calculations. 

For the M,($(i?e/8)) predictor, the coefficient of the quadratic term (ai) shows large 
variations with morphological sample. It varies from —0.21 ± 0.4 for the elliptical {15E) 
sample to 0.25 ± 1.06 for the spiral/lenticular {8B) sample to 0.43 ± 0.48 for the combined 
galaxy [23gal) sample The combined galaxy sample fit is logM, = (8.10 ± 0.12) + (1.71 ± 
0.24) log $(i?e/8)„ + (0.43 ± 0.48)[log<l>(i?e/8)„]2 (intrinsic uncertainty 0.36). As with the 
log-linear expression, this is not a strong predictor for galaxies as a whole because of the 
strongly differing elliptical and spiral/lenticular galaxy populations. 

For the M,(M(10i?e)) predictor, the coefficient of the quadratic term (ai) is always 
consistent with zero. The quadratic and linear (02) terms are larger for the spiral/lenticular 
sample than for the elliptical and combined galaxy samples (0.18 versus 0.08 for oi and 
1.0 versus 0.87-0.89 for 02). The combined galaxy sample fit is given by logM, = (8.02 ± 
0.09) + (0.87 ± 0.09) log M(10i?e)„) + (0.08 ± 0.08) [log M(10i?e)„)]2 (intrinsic uncertainty 
0.2 9). Thes e coeff icients are barely consistent, within the limits of the uncertainties, with 



the IWyithd (120061 ) result, although this may be, in part, because of the slightly different 
normalization of the mass. (Unlike in a log-linear fit, selecting a different normalization for 
the log-quadratic fits will result in different coefficients.) 



Although none of the predictive fits for any of the fitting samples showed a strong 
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indication for the necessity of a log-quadratic, as opposed to a log-linear, fit, neither do any 
of the predictors refute th at a quadratic does fit the data, and that the fits are compatible 
with the results quoted by IWyithd (120061 ) . For all of the predictive fits, the log-quadratic fits 
were statistically equivalent to the log-linear fits, using the statistical indicator adopted here, 
the Snedecor F-test. This i s a different statistical indicator than the Bayseian techniques 
employed by IWyithd (120061 ) . The M,(X) residuals versus X illustrated positive residuals 
for the three largest elliptical galaxies, but matching positive residuals are not seen for 
the small number of low mas s galaxies included in the our fitting sample. Two of the 
smallest dispersion galaxies in iTremaine et al.l (120021). NGC 474 2 and the Milky Way, are 
not used in our fitting samples, but are included by IWyithd (120061 ). It is possible that with a 
revised sample including more high and low mass galaxies that the need for a log-quadratic 
expression would become more evident using the statistical criteria employed in this paper. 
However, given the current sample, the simpler log-linear expressions are preferred over the 
log-quadratic expressions for all predictors. 



B.2. Multivariate Fits 

In addition to considering the log-quadratic expression, it is important to explore the 
M,(X, i?e) multivariate predictors for X = Eg,^{Re/8), and M(10i?e)- When this multi- 
variate expression was explored for X = a, it was found that the addition of was an 
improvement for the spiral/lenticular galaxies {8B), but not for the elliptical galaxy {15E) 
sample. Examining the residuals for the M,{X) residuals versus Re for the gravitational 
binding energy (Figure [9]), the gravitational potential (Figure [TTj) and the bulge mass (Fig- 
ure [T^ , there is no evidence for a residual- i?e correlation for the binding energy or bulge 
mass predictors. For the ^(-Re/S) predictor, there is a hint of a correlation for elliptical 
galaxies filled stars) when using the §i?-fit and the 23gal fit, but not when using the 

better-fitting 15E-fit. This is an indication of the lower slopes producing poor predictions 
of M, for elliptical galaxies, rather than a fundamental correlation. 

The tabulated multivariable M,{X,Re) fits for the gravitational binding energy, gravi- 
tational potential and bulge mass are given in Tables [191 EOl and [21] respectively. These 
tables enumerate the fits expressed as log(M,) = oi log (X„) + 02 log (i?en) + d, such that 
for each of the four morphologically based samples {15E, Cap8, 8B, and 23gal), the best-fit 
coefficients (a^, d) are given along with the requisite intrinsic uncertainty, the scatter, the 
maximum and minimum ai coefficient obtained by varying the weightcent selections, and 
the F-test ratio and its statistical significance (see Section [1]), both when compared with the 
log- linear M,(cr) predictor and when compared with the log- linear M,[X) predictor. 
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For the M,{Eg,Re) predictor, the single variable, log-linear M,{Eg) predictor still ap- 
pears to be the most reliable. For the elliptical samples, the M,{Eg, Re) predictor is not 
statistically significantly better than the M,{Eg) predictor, with a predicted fit of M, oc 
^o.88±o.i2^-o.62±o.26 (^5^), The spiral/lenticular sample {8B) suggests a relation of M. oc 
^o.4i±o.i7^o.59±o.76^ which is incompatible with the elliptical fit. The combined sample, how- 
ever, still suggests that M. oc E'^-^ (M. oc £;o.57±o.io^o.o2±o.24^ ^-^^^i little contribution from 
Re- Given the non-statistically-significant improvement for the elliptical galaxies when Re is 
added and the discrepant spiral/lenticular and elliptical predictors with the inclusion of Re, 
the single- variable M,{Eg) predictor will be taken to be superior to M,{Eg, Re) for galaxies 
in general. 

The M,($(i?e/8), Re) and M,($(i?e), Re) predictors show a slight, but not statistically 
significant, decrease in the residuals, when compared with M,($(i?e/8)) and M,($(i?e)) re- 
spectively. However, although the scatter is shghtly lower, the fits are still irreconcilably 
different for the elliptical (M. oc [$(i?e/8)]^-^^^"-2^i?°-22^°-^^) and the spiral/lenticular (M. oc 
[$(i?e/8)]°-^2±°-2^i?°-^^±°-^2) populations. The combined sample fit of M. oc [$(i?e/8)]^-^^^°-22^°-^^±°-^^ 
fits neither population particularly well. Thus, since there was no substantial improve- 
ment to the predictors, nor were the discrepant spiral/lenticular and elliptical predictors 
brought into closer functional agreement, the M,($(i?e/8)) predictor will be preferred over 
M.($(i?e/8),i?e)- However, the spiral/lenticular fit implies M. oc [^{Re/S)f-^'^R^e^^ ~ 
^0.82^0.11 ^ M°-^^. This, in turn, implies that a quantity akin to mass is a strong pre- 
dictor of M, for the spiral/lenticular galaxy population. 

Finally, the M,(M(10i?e), -Re) predictor shows a significant improvement (lower residu- 
als) over the M,(M(10i?e)) predictor for the eUiptical samples. The multivariate fit of M, oc 
[M(10i?e)]^-^^^°-^^-Rr^-^^^°-^^ for the 15E eUiptical sample suggests that a physical quantity 
similar to energy (M, oc M^-^'^R'^-^'^ ^ M'^R~^ oc £^) is a better predictor than pure mass for 
the elliptical galaxies. For the spiral/lenticular sample, M, oc [M(10i?e)]°-^^=^°-^^-Re^^'^°'^°' 
the contribution from Re is minimal and the scatter is (non-statistically-significantly) worse, 
suggesting that mass alone may be an an adequate predictor of M,. The "pure" mass- 
predictor for spiral/lenticular galaxies was also suggested by the M,{^{Re), Re) predictor. 
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Table 1. Surface Brightness Profile Sources 



Galaxy 


Ref 


Form 


Range 


H/G 


Color(Ref) 


Ficon 


(Ref) 


NGC 221 


1 


D 


0.02-1.1 


H 










2 


D 


0.85-116 


G 


ri=P(2)=' 


MN:e 


=P(2) 


NGC 224 


1 


D 


0.01-10.5 


H 










3 


D 


40.78-176.9 


G 


vg=-0.03-0.37gr;gr=P(3) 








4 


Fd 


130-5600 


G 








NGC 821 


5,6 


D 


0.2-327.9 


H,G 








NGC 1023 


7 


D 


0.2-8.7 


H 










4 


Fd 


6-284 


G 








NGC 1068 


8 


D 


0.7-11 


C 










4 


Fd 


12-280 


G 








NGC 2778 


5,6 


D 


0.02-348.5 


H,G 








NGC 2787 


9 


F 


-18 


H 


bv=l. 06(10) 


MN;€= 


=0.23(9) 




4 


Fd 


19-126 


G 


bv=1.06(10) 






NGC 3115 


7 


D 


0.2-3.9 


H 










4 


Fd 


5-360 


G 








NGC 3245 


11 


D 


0.01-3 


H 




MN:e 


=P(11) 




4 


Fd 


4-110 


G 


br=l. 39, bv=0. 87(12) 






NGC 3377 


5,6 


D 


0.02-367.1 


H,G 








NGC 3379 


13 


D 


0.02-14.6 


H 










14 


D 


0.4-4.8 


G 


vr=0.89(18) 








15 


D 


6-170 


G 


bv=0.96(10) 








16 


D 


8-255 


G 


bv=0. 96(10) 








17 


D 


8.3-159 


G 


bv=0. 96(10) 








18 


D 


8.3-191 


G 


bg=0.65bv(18),bv=0.96(10) 








19 


D 


8.3-347 


G 


bv=0.96(10) 


MN:c/ 


i=P(19) 


NGC 3384 


20 


D 


O.O2-I4.6 


H 










4 


Fd 


4-184 


C 








NGC 3608 


5,6 


D 


0.02-321.6 


H,G 








NGC 4258 




D/Fd 


0. 04-149. 4 


H,G 








NGC 4261 


14 


D 


1.3-4.8 


G 


vr=0.91bv(14),bv=0.99(10) 








23 


F 


0.1-10.25 


H 










15 


D 


10.6-140 


G 


vr=0.62(24),bv=1.0(10) 






NGC 4291 


5,6 


D 


0.02-327.9 


H,G 








NGC 4342 


25 


Fd 


0.01-5.51 


H 


vi=l. 3(25,26) 


MI:e=0. 


275 (25b) 


NGC 4459 


27 


D 


0.3-9.4 


G 










19 


D 


1.4-10.6 


G 




MNx/ 


i=P(19) 




4 


Fd 


10-128 


G 


bv=0.98(10,28) 






NGC 4473 


5,6 


D 


0.05-355.5 


H,G 








NGC 4486 


29 


D 


0.02-19.5 


H 


vi=1.73(29) 








27 


D 


1-983 


G 


bv=0. 96(10) 








19 


D 


1.4-84.7 


G 


bv=0. 96(10) 


MNx/ 


i=P(19) 




17 


D 


1.6-138 


G 


bv=0.96(10) 








15 


D 


4.6-163 


G 


bv=0.96(10) 






NGC 4564 


30 


F 


0.2-18 


H 


bv=0.89,vi=1.14,br=l. 44(31) 








32 


F 


15-119 


G 


bv=0.89,vi=l. 14(31) 






NGC 4596 


4 


Fd 


3-146 


G 


bv=0.96(10,28) 







Table 1 — Continued 



Galaxy 


Ref 


Form 


Range 


H/G 


Color{Ref) 


i?corr(Rcf) 


NGC 4649 


5,6 


D 


0.02-381.9 


H,G 






NGC 4697 


5,6 


D 


0.02-382.2 


H,G 






NGC 4742 


7 


D 


0.2-24 


H 


bv=0. 81(10) 






33 


Fd 


2-80 


G 






NGC 5845 


5,6 


D 


0.02-342.3 


H,G 






NGC 6251 


34 


D 


0.03-1.1 


H 


ubpoc = l-Subj + 0.4;ub=P{34),bv=1.0(10b) 


MN:e=P,0.112(35) 




35 


D 


1.5-6.1 


G 


vr=0.91bv{14);bv=1.0(10b) 


MN:e=P(35) 




35 


D 


2.2-18.4 


G 


brg = 0.34 + 1.57to(36);bv=1.0(10b) 


MN:e=P(35) 


NGC 7052 


14 


D 


0.4-2.2 


G 


vr=0.91bv(14);bv=1.0,vi=1.4{37) 






37 


D 


2.5-15.8 


H 






55 


38 


D 


6.7-31 


G 


brg = 0.34 -1- 1.57feD(36);bv=1.0,vi=1.4{37) 


MN:€=P(38) 




39 


D 


28-119 


G 


vi=1.4{37) 




NGC 7457 


20 


D 


0.03-5.5 


H 








40 


Fd 


6-40 


G 


br=1.287,rk=2.403(41),bv=0.89{10) 


MN:e = 0.35(40) 


IC 1459 


42 


D 


0.05-10.6 


H 


vi=A{43) 


MN:e=P,0.263(42) 




44 


D 


1.5-77 


G 


br=P{44);bv=0.98{10);vi=A(43) 


MN:£=P{44) 



^For conversion to V for SB pro file comparisons the follow ing conversions were used; N221 (vi=1.133 l lTonrv et aLllioa 
and NGC 3245 (br=1.39,bv=0.87 l lMichard fc PoulainllSoOol) '! 



Note. — Description of columns: (1) galaxy name; (2) SB profile reference; (3) the form of the data: F (fit such as 
Nuker-law, de Vaucouleurs or Sersic profile), Fd (fit for the bulge portion of a bulge-disk decomposition), or D (original data 
points); (4) radial range (in ") along major axis; (5) ground (G) or HST (H) based data; (6) color correction(s) applied to 
original SB profile such that a color correction of X-Y will be given as 'xy'. For those cases where a single color correction 
was used for multiple radii, that is given. For those galaxies where the color correction was a function of radius, the value is 
listed only as P here, and for those where there were several values within specified aperture radii the value is listed as A, 
with the source of the correction; (7) radial correction(s) applied to correct the original SB profile to the major axis listed 
such that MN implies a correction from the mean axis and MI implies a correction from the minor axis, with the ellipticity 
(e) or axial ratio (q) used in the conversion given such that if a single value was used it is stated, or if the ellipticity as a 
function of radius was used, this is given as P, with the source of the correction. The galaxies not belonging to the 23gal 
sample, i.e. those galaxies not used in the primary fitting functions, are denoted in italics. 



References. — (D iLauer et al.l l ll998h . f 2) IPeletieij lll993l ). (3)'Hoessel & MelnicW (ToSo') - "positions with respect to th e cen- 
ter ha ve been multiplied by a factor of 1.035" as given by [Waltcrbps & Kcnnicutt (1987), (4) Baggctt, Baggctt, & An dersonI 
lll998l'). (5) Full profile obtained directl y fro m K. Gebh a rdt, (6) Ictebhardt et al.l ||2003|). f 7) iLauer et al.l lll995l 'l. (8) 
ISchild. Tresch-Fienbere. fc Huchral l ll985h . (9) IPeng et al.l (120021 '). ( 10) Ide Vaucouleurs et al.l (Il991), (10b) typical value 



de Vaucouleurs et al 



■ ( ll991^ for other sample ga laxies , (11) take n fro m Figure 2 of . . .. . , , ^ 

iMichard fc PoulainI boodt) . fl S*) ICebhardt et al l toOOs ), (MllLaueij (ll985^ . (ISllPeletier et al.l (Il990h . ('16') lBursteinr (Il979l '). 
(17) iDavis et al.l ( Il985ir ^fl8') iKormendvl (Il977f). (19) llVIichardl (Il985h. (20) [ Lauer et al.l (I2005I') . (21) Full profile obtained 
dire ctly from C. Siopis, (2 2) 'Siopis et al.l (|2006|), (2 3) 'Ferrarese et al.! ('199I) with 0.5 mag added to m atch other data, 
(24) iFerrarese, Ford, fc Jaffc ( 1996). (25~ IScorza fc van don Bosch (1998). (25b) IScorza fc van den BoschI jlgO S) Figure 1, 



Barth et al.] ||200]J), (12 ) 



(34) 
(35 



26) ICretton fc van den BoschI (Il999[). (27)ICaon. CapaccioU, fc Rampazzol lll990l). (28) Sarzi et al.l ( I200J). (29)lLauer et alj 



1992a^. (30)lRest et al.l 1I2OOII ). (31) lMichardl ( l200(t ). (32) ICaon. CapaccioU. fc D'Onofrid ( Il993h . (33) 1vigroux et al.l (Il988l ). 



Grane et al.l lll993l) Figur e 23 assuming (U — B) ppc — 1.8((7 — i?) j+0.4 may be ap plied to all radii and that Upoc ~ U. 



Young et al.l ( Il979l). (36) Ide Juan. Colina, 



_ _ , , , „ Perez- FournonI lll994l); iKenj 1^^85 1). (37) [vai^de^Mare^^^ van den BoschI 

lIlQQsl) Figure 3. (38) Ide Juan. Colina. fc Perez-FournonI (Il994l ). (39) lGonzalez-Serrano fc Perez-FournonI (ll992^ F i gure 6 , (40) 
iKhosroshahi. Wadadekar. fc Kembhavil (|200(]|) with 0.79 mag added to match other data. (41) |Peletier fc Balcellj (Il997l) , (42) 
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iForbes. Franx. fc niineworthl lIlQQSl '). (43) ICaroUo et all (Il997h . (44) iFranx. Dlingworth. fc HeckmanI l ll989h 



Table 2. Host Galaxy Parameters Used in Modeling 



Galaxy 




Af.{low,high)'^ 


Ref 






Ref 




Band 


q 


Ref 


e 


Ref 




Ref 


Sample 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 


Mr'n 991 




z.yiiiDi^z, 0,0.0 ) 


1 
1 


fl SI 
U.Ol 


1 


1 
1 


1 '^7 
1.0 ( 


T 
1 


n 7'\ 

U. 1 


1 
1 


70 
1 U 


1 
1 


n 078 

u.u / 


9 
Z 


1 f^TT Por^S 

ioi!j,'^apo 






4. oil; / l^Z.U,o.O ) 


Q /I c; c 
0,4, 0,D 


n 7fi 
U. / D 


O.U 


D 




V 


U.ozO 


^ 7 
0, / 


77 


Q fi 
0,0,0 


1 (\A 
1.U4 




y 


OJj 


Tsirr* S91 






1 n 


94 1 
Z4. 1 


7 R 


1 n 


zl ni 
4.U1 


V 


n fi9 
u.uz 


1 1 


yu 


1 n 

lU 


(\ 7(\ 


1 n 

lU 


15E,Cap8 


INvj*^ lUZo 


orsi^rsju- 


4.411; I i^o.y ,4.0 ) 


1 9 
iZ 


i i.4 


/I 07 
4.y ( 


1 9 

iZ 




V 


n 79 
u. /z 


1 9 

iZ 


on 
yu 


1 9 
iz 


n A Q/1 



y 


OO 


iVLrC' lUOo 


( JT jo/±(rsjo 


1 . 1 [1 .U,o.U/ 


0,10 


10. U 


l.Zl 


y / 

-^4 




V 


U.o fO 


1 


00. 


1 


1 1 Q 

1 . lo 


n 

y 




iV LrO Z 1 lO 




1 [U.O,Z.Z/ 


1 u 


ZZ. bf 


0. U 


1 U 




T/ 

V 


u. /o 


1 1 
1 1 


on 


1 u 


z. uo 






ivrpr^ 97S7 
IN V_i V_ Z ( c5 1 




4.1il;i l^O.D,4.0j 


1 7 


7 ^ 
1 .0 


91 
0. Zl 


1 s 




■R 
Jj 


n 77 


1 Q 

ly 


OU 


1 7 
1 1 


U.DoU 


Q 

y 


8R 
OJj 


INVjr^^ O±±0 


OU- 


1 nTTQl'n A (W 
±.Ull;y l^U.4,Z.U ) 


p. on 91 

D,ZU,Z1 


Q 7 

y. / 


D.y 


D 




V 


U.D 


91 
Zl 


oD 


91 99 
Zl ,ZZ 


9 87 
Z.o / 


Q 

y 


OJj 






Z.lll;Ol^i.D,Z.D ) 


9Q 
Zo 


on Q 
zu.y 


7A 
0. ( 4 


9Q 
ZO 




rt 


n fit; 
U.oO 


O'XV^ 
Zv5D 


Do 


9*5 


i.OO 



y 


SR 
OJj 


IN V_i V._ OO / 1 


IIjO-D 


±.Ull;c5l^U.y,l.y ) 


1 n 


119 
1 ±.z 


9 Q 

z.y 


1 n 


9 RQ 

z.oy 


V 


u.ou 


1 1 
1 1 


on 
yu 


1 n 

lU 


9 1^7 
Z.O / 


1 n 

lU 


ioii;,'^apo 


Nir^r^ Q07Q 




i.Uil;Ol^U.Di,Z.U ) 


9/1 
Z4 


1 n fi 


A 7fi 
4. i D 


9/1 
Z4 


/I 


V 


n on 
U.yu 


9C^ 
ZO 


on 
yu 


OA 
Z4 


•3 9n 
o.zU 


9(^ 
ZD 


10ll;,L.'apO 


iVLrC' 0004 


OlJ[SJU- 


1 . orj 1 [1 .4, J- ■ 1 ) 


1 u 


11.0 


Z. 


1 u 




T/ 
V 


U. OU 


1 1 


on 


1 U 


U. loO 


n 

y 




IN^O oDUo 


T?9 


i.yii;oi^i.v5,z.y J 


lU 


99 

zz.y 


Q 7 


1 n 


4. DO 


V 


n Q9 
U.oz 


1 1 


on 
yu 


1 n 
iU 


A 0^ 
4.Z0 


1 n 

lU 


1011;,L.'apO 




Or\.D[S J DC 


0. (74 -C/ 1 [ <■. ifi J . Z4 J 


Z 1 


1 . Zo 


q q 


/ 






u. 00 


^ / 


IZ 


z 1 


14 - 00 


0-7 




■pvT/^r^ /19fi1 


T?9 


^ OT?'>i( A A « 0\ 
D.Zll;ol^4.1,0.Z ) 


9fi 
ZO 


oi.D 


A 7^ 
4. / 


9fi 

zo 




V 


n 7^ 
U. / 


9Q 90 

zo,zy 


K/1 

D4 


9Q 
ZO 


nn 
o.UU 


9R 

ZD 


ivrpp 49Q1 

iNvj*^ ^zyi 




0. iIl;ol^U.c5,o.y ) 


1 n 


0(<. 

ZD.Z 


0.0 


1 n 




V 


n 7(^ 


1 1 
± 1 


on 
yu 


1 n 

lU 


9 19 

z. iz 


1 n 

lU 


lOH; 




o(/- 


Q 1 77Q /<? n / S) 
o.lrjo[ Z.U,4-oJ 


(/ 


1 K Q 
10.0 


O.Z 


oU 




J 
1 


U. lO 


Q 1 
01 


iJU 




U. ZUU 


z 




iNvx*^ 4^oy 






1 7 


ID. i 


D.o 


•^9 
oZ 


A 99 
4.ZZ 


13 


U.oO 


00 


A^ 

41 


1 7 
1 1 


1 '^n 

l.OU 


Q 

y 


SR 
OJj 


NGC 4473 


E5 


1.1E8(0.31,1.5) 


6 


15.7 


6.0 


10 


3.66 


V 


0.62 


11 


72 


10 


2.34 


10 


15E,Cap8 


NGC 4486 


E+0-1 


3.0E9(2.0,4.0) 


6,34,35 


16.1 


4.0 


6 


8.27 


V 


0.915 


36 


42 


35 


8.59 


37 


15E,Cap8 


NGC 4564 


E6 


5.6E7(4.8,5.9) 


10 


15.0 


2.0 


10 




/ 


0.70 


11 


90 


10 


1.84 


10 




NGC 4596 


SB{r)0+ 


7.8E7(4.5,11.6) 


17 


16.8 


1.2 


38 




B 


0.85 


39 


42 


17 


1.42 


9 


8B 


NGC 4649 


E2 


2.0E9(1.4,2.4) 


10 


16.8 


8.5 


10 




V 


0.88 


11 


90 


10 


6.34 


10 


15E 


NGC 4697 


E6 


1.7E8(1.6,1.9) 


10 


11.7 


4.7 


10 




V 


0.6 


11 


90 


10 


5.49 


10 


15E 


NGC 4742 


SO 


1.4E7(0.9,1.8) 


6,40 


15.5 


2.67 


26 




B 


0.6 


41 


58 




0.339 


^5 




NGC 5845 


E3 


2.4E8(1. 0,2.8) 


10 


25.9 


5.5 


10 


4.54 


V 


0.71 


11 


90 


10 


.605 


10 


15E,Cap8 


NGC 6251 


E2 


5.3E8(3.5,7.0) 


44 


93.0 


9.7 


44 




V 


0.85 


45 


76 


44 


6.81 


2(&45) 


15E 


NGC 7052 


E4 


3.3E8(2.0,5.6) 


46 


58.7 


6.3 


6,46 




I 


0.55 


47 


70 


46 


6.69 


48 


15E 


NGC 7457 


SA(rs)0- 


3.5E6(2.1,4.6) 


10 


13.2 


3.2 


10 


2.13 


V 


0.7 


11 


90 


10 


0.972 


2 


8B 


IC 1459 


E3 


2.5E9(1.4,3.6) 


49 


29.2 


3.2 


49 




I 


0.725 


50 


90 


49 


6.75 


49 


15E 



'"Galaxy type s are fro m NED for all galaxi es except for NGC 2778, NGC 4742 NGC 5845, NGC 62 5 1 and NGC 7052. The t ype given here for those galaxies is from 
iGebhardt et al.i (2003'): Tremaine et al.l i l2002l) for NGC 2778 and NGC 5845, from lTremaine et al.l ll2002h : iFerrarese fc FordI lll999l) for NGC 6251, a.nd fr om lTVemaine et all 
l l2002h : lvan der Marel fc van den BoschI lll99a) for NGC 7052. NGC 4742 is classed as SO in lLauer et all lll995l ') but as E4 in NED lTremaine et al.l ll2002h . 



^Galaxy classific ation in the literature differ s from NED as follows: NGC 821 - classifie d as E4 dOebhardt et al.ll2003l: I ^^m^me et al. 200 ^), NG C 3115 classed as an E 
l lLauer et al.lll99d). NGC 4291 classed as E2 I Gebliardt et al.ll2003l : iTremaine et al.ll2002l') NGC 4342 noted as E7 |Scorza fc van den Boschlll998h , NGC 4564 classed as 
E3 iGebhardt et alj 20031: iTremaine et alJbOOa) . NGC 4649 classed as El llGebhardt et aUl2003l: iTremaine et al.ll2002l) and NGC 4697 classed as E4 llCebhardt et al.ll2003l: 
iTremaine et al.ll2002l') " 

"^The values in ITremaine et aL I 1 I2OO2I) for the black hole masses (and limits) differ from those here for the following galaxies. The value in ITremaine et al J l l2002h : NGC 
221 (2.5x10'' (2.0,3.0)), NGC 1023 (upper limit of 4.9), NGC 3379 (5.0x10'' and 1.6x10^ as limits), NGC 4258 (3.9xl0''(3.8,4.0) from lHerrnstein et al.l lll999l) ). NGC 4342 
(3.0x10® (2.0,4.7)), NGC 4596 (upper limit of 12.0), IC 1459 (limits (2.1,3.0)). 

''The distance for each of the galaxies is taken directly from ITremaine et al. I ll2002h . llTremaine et al. I l l2002h used SBF distances l lTonrv et al.ll200lh . when available, and, 
when not available, from the redshift assuming a Hubble constant of 80 kms~^ Mpc~^ .) 

''The color terms used for conversion to the default wavelength-band arc as follows with V-I(vi) colors from lTonrv et al.l 1 I2OOII) and B-V(bv) color from lSarzi et al.l ^0^: 
NGC 821 vi=1.196, NGC 3377 vi=1.114, NGC 3379 vi=1.193, NG C 3608 vi=1.156, NGC 4459 vi=1.187, bv=0.98, NGC 4473 vi=1.158, NGC 4486 vi=1.244, NGC 5845 
vi=1.124, & NGC 7457 vi=1.104. The solar colors were taken from iLivinestonI ||2000| ) such that b?;0=O.65O, ufe0=O.195, ui;0=O.845, ^©=0.54, ■!;i0=O.88, ■!;A:q=1.49. 

teender et al.l 1I2OO5I ) give a higher SMBH mass of 1.4x10* Mq (1.1,2.3), based on HST STIS spectroscopy and a revised model for the nucleus of NGC 224. The implications 
of using this alternate mass are discussed in Section |4] 



Note. — Description of colum ns: (1) galaxy name; (2) galaxy type from NED; (3) black hole m ass (in M(Ti) from (4) re ference, rescaled to (5) distance (in Mpc); (6) 
best-pre- ICappellari et al] 1 I2OO6I ) mass-to-light ratio, Tjpc (in Mq/Lq), from (7) reference, and (8) ICappellari et al.l ll2006h mass-to-light ratio, Tcap, (in Mq/Lq), both 
in (9) wavelength-band, at the distance given in col. 5; (10) observed axis ratio from (11) reference; (12) inclination (in °) from (13) reference; (14) bulge effective radius 
along the major axis (in kpc) from (15) reference; (16) sample membership denoted as 8B, 15E and/or Cap8. All members of these sub-samples also belong to the 23gal 
sample. The galaxies not belonging to the 23gal sample, i.e. those galaxies not used in the primary fitting functions, are denoted in italics. 
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References. — (DlVerolme et al.l ll2002h . (2) crude r^/* law fit to composite SB profile, (3) 'B acon et al.l teOOlV (4) lKormendv fc Bender 199E ). (5) l^emaini |[l99l) , (6) 



ITremaine et al.l ll2002h . f7) lKenl] lll989l). (8) IWalterbos fc Kennicutt 
tion with K. Gebhardt, (12 ) Bower et ahl l l200l[). Cl3) 



Greenhill et al 



fl987), (9) Baggett. Baggett. fc Anderson 119981 ), (10) ICebhardt et al ll20q3l). (11) private communica. - 
1199a ). (14)jTkkajnija ir^ofa^l l200Cl) Figure 4 inner regions, (15) average of lde Va u£ouleurset aL 1 199ll ) 



(0.85) and lMarcum_ et"al| 1I2OOII) (0.90). (16) ISanchez-Portal et al.l ^20001) . ( 17) ISarzi et al.l 1I2OOI \, ( 18) Calculated using (B) L^^ige = 3.25xlO^LfT^ llSarzi et al.ll200ll) and 



iBinnev fc Tremainj 1119871 ) eq. 4-80b (M ~ < v'^ > r^/ (0.4G)) with = 184.9fcmg~^ jSarzi et all200lli to obtain = 3.38. A second calculation wherein M = v^^^R^jG 



Barth et al.l ||2001 



Figure 2 of 

lll985h :|de Vaucouleurs, 



), (24) 



Gebhardt et al. (2000a) 



ICaon^ C apaccioli. fc D 'Onofrio 1 1994):IScorz a 



Macchetto et 



al.l 1I1997I) . 



(^'m gj; = V sm Q') and using Q = 50° jSarzi et al.ll20pl ) an d Vmax = ISOfcmg ^ f rom Figure 4 of Sarzi et al. (2001), result ed in = 3.03. T hese were averaged to 3.21., 
(19) IPeng eTa l. (2002j^ (20) iKormendv et al. l l996h . (21) i Kormendv fc Richstonc (19J2), ( 22) Capac cioli,^ H eld, fc Nic toi lll987 | ). (23) [Sarth e t al. (2 001), (23b) Based on 



(25) central value o f iFaber elTall (Il989l): iBender. Dobereiner. fc MoUen hoff (1988j): iBaconTMonnet. fc Simieiil 



Capacciol: (Il979l) : Burstein (,197£ ): Peletier et al.l llT99ol): van der Marel. Binney fc Davied (ll99C 
SimienI 1 198n . (27) ISiopis et al l ll200e ), (28) iFerrarese. Ford, fc Jafi'j 1 1996 ). (29) Ferrarese et al.l 



van den BoschI lll998l) . (32) lFaber fc JacksonI 197£ ). (33)(Michar dl (^igSSl) 
(35) lHarms et al. (1994), (36) central value of lFaber et al.l (Il989l) : IPeletier et al.l 61990) 
(38) Calculated, as in (18), using Lfcuige = 9.61a:10®Z/Q, ere = 136.3fc mg~^, iimaa; = 80fcms~', and d = 42° (Sarzi et al 



van der Mareli (il99ll) values which are 0.89 -0.91, 



19941). (30)ICretton fc van den Boscl] ^99^, (31) 



Cagn^Cagaccioli^^^^mgazz d ill 990^ profiles, 
Ba con, Monnet. fc SimienI l ll985^ (37) |Peletier et al.l 
20011 ). and obtaining T g = 1.47 in the first 



calculation and = 0.91 for the second, with an average value of 1.19., (39) IWozniak fc Pierce! l ll99d) SB pr ofile, (40) Trcmaine et al.l(l2002[) (whic h referen ces the unavail - 
able M.E. Kaiser et al. 2002, in preparation) (41) based on lLauer et al.l lll995h SB profile, (42) using method in lYasuda. Okamura. fc Fukugita lll995h based on lHubbl3 (Il926l ) 

(44) 'Ferrarese fc Ford] 
Perez-Fournon. 1199^) . 



{i = cos-^^J{ql^ - 0.22)/(l - .22)) which here would result in 58°with the'de Vaucouleurs et al.' ("19911) value of 925 = 0.56, (43)IVigroux et al.' ('l988^ 
l ll999h. (45)lBender. Dobereiner. "fc Mollcnho ff (1988). (46) Ivan der Ma rel fc van den Bos ch (1998), (47)lFaber et al.lll989 | ). (48) Gonzalez-Se rrano fc 
(49) ICappellari et al.l (I2OO2I) . (50) average of IPaber et al.l lll989h ESQ. FW and iFranx.lUingworth. fc HeckmanI l ll989l) : IVerdoes Kleiin et al.l (l200ol) ~ 



Table 3. Galaxy Parameters Used in the Fits 



Galaxy 


cr 




M{10Re) 


M{Re) 


M(Wkpc) M{kpc) 


*(i?e) 


<I>(R,;/8) $(20pc) 


${100pc) 


Eg 


1 


2 


3 


4 


5 


6 7 


8 


9 10 


11 


12 



NGC 0221 
NGC 0224 
NGC 0821 
NGC 1023 
NGC 1068 
NGC 2778 
NGC 2787 
NGC 3115 
NGC 3245 
NGC 3377 
NGC 3379 
NGC 3384 
NGC 3608 
NGC 4258 
NGC 4261 
NGC 4291 
NGC 4342 
NGC 4459 
NGC 4473 
NGC 4486 
NGC 4564 
NGC 4596 
NGC 4649 
NGC 4697 
NGC 4742 
NGC 5845 
NGC 6251 
NGC 7052 
NGC 7457 
IC 1459 



75 

160 

209 

205 

151 

175 

140=* 

230 

205 

145 

206 

143 

182 

105 

315 

242 

225 

186 

190 

375 

162 

152 

385 

177 

90 
234 
290 
266 

67 
340 



[-02 
[-02 
hOl 
hOl 



2.52E-(-03 
4.12E-f02 
2.95E-f01 
1.70E+03 
2.11E+03 
7.22E+01 
3.87E-(-02 
3.18E-i 
3.99E-i 
9.49E^ 
9.12E^ 
7.39E+02 
5.00E-(-01 
2.13E+01 
6.39E+01 
1.76E+02 
2.43E+03 
4.33E-f02 
1.55E-f02 
6.36E+01 
2.24E+O2 
3.28E-(-02 
1.04E+02 
5.75E+01 
3.06E+03 
1.12E+03 
4.73E+01 
7.00E-f01 
1.02E+02 
1.03E+02 



3.48E+08 
3.44E-flO 
6.02E-(-10 
2.18E+10 
3.96E+10 
2.12E+10 



[-09 
-11 
-10 
-10 
-10 



3.59E^ 
2.29E-i 
6.29E-i 
1.86E^ 
7.54E^ 
I.43E+IO 
7.80E-(-10 
1.83E+11 
3.15E-fll 
5.61E+10 
1.33E+10 
2.64E-flO 
4.02E-flO 
7.81E+11 
2.33E+10 
7.78E-(-09 
6.62E+11 
1.25E+11 
5.13E+09 



2.26E- 
3.64E- 
2.85E- 
2.95E- 
3.85E- 



10 
-11 
■11 
-09 
11 



1.37E-I-08 
1.24E-I-10 
3.05E-I-10 
7.43E-I-09 
1.50E+10 
1.06E+10 
1.05E-I-09 
9.55E-(-10 
2.73E-(-10 
9.20E+09 
3.77E-I-10 
5.76E+09 
3.82E+10 
7.24E+IO 
1.29E-I-11 
2.61E-I-10 
5.75E+09 
9.51E-(-09 
1.69E-I-10 
2.55E-I-11 
9.70E+09 
3.21E-I-09 
2.81E-I-11 
6.19E-I-10 
1.60E+09 
9.51E-I-09 
1.52E-I-11 
8.84E-I-10 
1.37E-I-09 
1.64E-I-11 



3.60E-I-08 
3.46E-(-10 
4.13E-I-10 
2.32E-I-10 
4.O2E+IO 
2.00E+10 
3.88E-I-09 
1.95E-(-ll 
5.99E-fl0 
1.74E-I-10 
6.45E-I-10 
1.50E+10 
5.93E+10 
6.61E+10 
1.63E-(-ll 
5.14E-I-10 
I.4IE+IO 
2.63E-I-10 
3.65E-I-10 
3.04E-I-11 
2.25E+10 
7.67E-(-09 
3.95E-I-11 
9.26E-I-10 
5.73E+09 
2.33E+10 
2.07E-I-11 
1.61E-I-11 
2.99E-I-09 
2.23E-I-11 



3.55E-f08 
1.31E+10 
9.25E+09 
1.49E+10 
1.65E+10 
6.86E+09 



-09 
-10 
-10 
H09 
-10 



1.96E-I 
5.02EH 
2.12EH 
5.94EH 
1.73EH 
7.92E+09 
1.50E-flO 
6.82E+09 
2.65E+10 
1.72E+10 
1.22E+10 
9.48E-f09 
1.13E+10 
3.16E-flO 
6.61E+09 
2.97E+09 
5.37E+10 
1.71E+10 
4.32E+09 
1.48E+10 
2.32E+10 
2.17E+10 
1.57E+09 
4.16E-I-10 



1.53E-23 
1.03E-22 
3.26E-23 
1.57E-22 
1.18E-22 
3.75E-23 
1.82E-23 
2.75E-22 
1.31E-22 
2.79E-23 
8.11E-23 
6.21E-23 
6.26E-23 
4.19E-23 
1.28E-22 
9.00E-23 
2.23E-22 
6.77E-23 
5.98E-23 
2.62E-22 
4.43E-23 
1.87E-23 
3.37E-22 
8.36E-23 
5.51E-23 
1.28E-22 
1.72E-22 
1.43E-22 
1.07E-23 
1.87E-22 



3.75E-23 
2.05E-22 
1.05E-22 
3.29E-22 
2.68E-22 
1.06E-22 
3.94E-23 
7.32E-22 
3.43E-22 
8.64E-23 
2.36E-22 
1.53E-22 
1.89E-22 
1.08E-22 
3.63E-22 
2.59E-22 
5.44E-22 
1.53E-22 
1.62E-22 
6.10E-22 
1.08E-22 
4.91E-23 
8.42E-22 
2.39E-22 
1.33E-22 
2.89E-22 
4.13E-22 
3.46E-22 
2.71E-23 
5.32E-22 



2.72E-23 
2.38E-22 
2.39E-22 
3.74E-22 
3.26E-22 
2.12E-22 
4.90E-23 
1.16E-21 
5.23E-22 
1.87E-22 
3.79E-22 
1.87E-22 
3.62E-22 
2.30E-22 
5.41E-22 
3.84E-22 
5.55E-22 
1.83E-22 
2.14E-22 
7.95E-22 
1.72E-22 
6.08E-23 
l.llE-21 
3.95E-22 
I.4IE-22 
3.32E-22 
6.69E-22 
4.50E-22 
3.49E-23 
8.91E-22 



1.14E-23 
2.09E-22 
1.87E-22 
2.74E-22 
2.74E-22 
I.4IE-22 
3.41E-23 
9.36E-22 
4.00E-22 
1.18E-22 
3.32E-22 
I.43E-22 
2.99E-22 
1.90E-22 
5.19E-22 
3.26E-22 
3.09E-22 
1.60E-22 
1.90E-22 
7.77E-22 
1.28E-22 
5.31E-23 
1.08E-21 
3.42E-22 
8.00E-23 
2.63E-22 
5.73E-22 
4.36E-22 
2.72E-23 
8.21E-22 



2.96E-15 
1.60E-12 
1.60E-12 
1.59E-12 
2.4OE-I2 
5.67E-13 
3.11E-14 
3.79E-11 
4.75E-12 
4.30E-13 
4.10E-12 
4.93E-13 
3.45E-12 
4.44E-I2 
2.46E-11 
3.26E-12 
1.63E-12 
9.04E-13 
1.47E-12 
l.OlE-10 
5.85E-13 
8.41E-14 
1.21E-10 
7.17E-12 
1.60E-13 
1.53E-12 
3.31E-11 
2.14E-11 
1.86E-14 
4.50E-11 



Normalization 200. 



200. 



4.0E10 



1.6E10 



3.6E10 



I.IEIO 



8.0E-23 



2.0E-22 



3.0E-22 



2.4E-22 



2.0E-12 



''The value of 140 kms ^ for cr, the rms stellar velocity dispersio n within a slit aperture of len gth 2i?e, is taken from iTremaine et al 
that sources in the literature quote a much higher central value, e.g.. lBarth. Ho. fc Sargent! ||2002| ) claim 202 fcms-l. lNeistein et al.l l ll999h 



||2002| '). It is noted 
claim 205 kms~^ , 



ISarzi et al. claim 210 kms ^, and lErwin et al" claim 257 kms ^. This discrep ancy between t he a v alue in | Tremaine et alJ | |2()0^ and the 

central value may stem from the central peak in the velocity dispersion profile, as illustrated in lBertola et al.l d 19951 ) . [Srwin e^dTil2003l ) and Neistein et al.l 
lll999l') . The implications of using a higher velocity dispersion, 200 kms ^, for NGC 2787 are discussed in Section HI 



Note. — The values of the bulge parameters used in searching for an alternative predictor of M, to the M, (cr) relationship are presented here for the (1) 
full sample of 30 galaxi es; not all 30 galaxies a re used in the fitting as discussed in the following se ctions. The (2) effective stellar bulge velocity dispersion 
in kms~^ is taken from lTremaine et al. 1 1 I2OO2I) for all galaxies except NGC 4258 llSiopis et al.ll20o3 ). The remaining parameters are in units of pc, Lq, Mq 
and seconds, or appropriate combinations thereof, with the remaining columns as follows; (3) V-band Ie\ modeled mass enclosed by sphere of radius (4) 
10-Re, (5) -Re, (6) 10 kpc, (7) Ikpc; modeled gravitational potential at (8) -Re, (9) Re/S, (10) 20 pc, (11) 100 pc; and (12) modeled gravitational binding 
energy. The bottom row contains the normalization applied to the data prior to computing the -M.-predictor fits. This normalization value is based on the 
mean or median (in log-space), or a compromise between the two values. The remaining normalizations not included in the table are: .Re (1600), T^^^^/^-R^ 
(3.2x10^''), I^Rl (1.8x10"), /e-Re (5.5x10^), and '^Ipc^l^l (3.98x10^5). The galaxies not belonging to the 23gal sample, i.e. those galaxies not used in 
the primary fitting functions, are denoted in italics. 
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Table 4. 


Quick Reference Table of Analysis Abbreviations 


Abbrv. 


Cat. 


Definition 


basic eq 


E 


Y = aiXi + d; Y = logM.; Xi = logx^; Xi=host gal. param. 


Isqfit* 


A 


multivariate least-squared fitting algorithm; utilizes ey 


fitexy 


A 


linear least sauares fittine alECorithmi utilizes e„ and e-r fPress et al. 1992") 


robust 


A 


unweighted robust (absolute deviation minimization) fittine algorithm (Press et al. (1992) medfit) 




WC 


uncertainty in Y 


fx 


WC 


uncertainty in Xi (default = 0) 


Avgerr* 


WC 


ty is symmetrized by averaging upper and lower Y-uncertainties 


Recent 


WC 


Y recentered between maximum and minimum Y (based on Y-uncertainties) 




WC 


"intrinsic" uncertainty in Y selected such that = 1 for fit 




WC 


"observed" uncertainty in Y 


OBS 


WC 


^ y — 


INT 


WC 




OBS+INT* 


WC 




weightcent 


WC 


Selection of Y-centering {Avgerr /Recent) and Y-weighting (OBS, OBS+INT, INT) 


N 


E 


number of host galaxy 'data points' in fit 


DOF 


E 


degrees of freedom in fit 




E 


number of 'new' data sets in bootstrap with Xr — 1-0 


23gal 


S 


final galaxy sample of 23 elliptical, spiral & lenticular bulges 


15E 


S 


15 elliptical galaxies in 23gal sample 


8B 


S 


8 spiral/lenticular bulges in 23gal sample 


Cap8 


S 


8 galaxy subset of 15E sample with Ycap available 




Q 


'^/y = DOF-^Y.^ {y^ - (Ej ^jXj, -f dif 


NSS 


Q 


difference between fits not statistically significant (< 75%) 


75SS 


Q 


difference between fits significant at 75-90% level 


90SS 


Q 


difference between fits significant at 90-95% level 


95SS 


Q 


difference between fits significant at 95-99% level 


99SS 


Q 


difference between fits significant at > 99% level 


'minimal improvement' Q 


scatter lower than reference (e.g. M,{a)), but not at 75% probability level 


strength 


Q 


a stronger fit has lower scatter 


stable 


Q 


a stable fit shows little variation with fitting algorithm, weightcent, and/or morphology 


best fit 


Q 


the most stable and strongest fit 



Note. — This table presents a reference to the abbreviations and codes (column 1) used throughout Sections |4} |5] The second column 
lists the general category of the abbreviated quantity: A=fitting algorithm, E= equation or variables therein, Q = quality of fit, S = 
fitting sample and WC=weighting/Y-centering. The final column gives a brief definition for the abbreviation. Starred parameters denote 
the default selections. 
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Table 5. Fitting Parameters: logM, = a log — + d 



Sample 


Algorithm 


a 


d 












15E 


Isqfit 


3.86±0.3S 


S.lSitO.OS 


0.23 


0.000 


0.28 


3.S9±0.41{AI) 


3.63±0.36(RO) 


23gar 


Isqfit 


3.75±0.32 


S.16±0.06 


0.24 


0.000 


0.30 


3.79±0.34{AI) 


3.71±0.32(RB) 


8B=- 


Isqfit 


3.27±0.77 


S.09±0.14 


0.30 


0.000 


0.38 


3.64±0.80(RI) 


3.27±0.77{AB) 


Caps 


Isqfit 


3.92±0.65 


S.16±0.14 


0.32 


0.000 


0.36 


4.00±0.49(AO) 


3.79±0.73(RI) 


15E 


fitexy 


3.90±0.39 


S.lSitO.OS 


0.22 


0.021 


0.28 


3.94±0.42{AI) 


3.S0±0.41(RB) 


23gal 


fitexy 


3.79±0.32 


S.16±0.06 


0.23 


0.021 


0.30 


4.01±0.1S(AO) 


3.75±0.33(RB) 


SB 


fitexy 


3.31±0.79 


S.09±0.14 


0.29 


0.021 


0.38 


3.69±0.81(RI) 


3.31±0.79(AB) 


Caps 


fitexy 


3.96±0.66 


S.16±0.14 


0.31 


0.021 


0.36 


4.07±0.2S(AO) 


3.S2±0.74(RI) 


15E 


robust 


3.99 


S.16 


0.00 


0.00 


0.28 




3.82(R) 


23gal'' 


robust 


3.70 


S.19 


0.00 


0.00 


0.30 




3.46(R) 


SB'' 


robust 


3.34 


S.13 


0.00 


0.00 


0.37 


3.48(R) 




Caps 


robust 


3.91 


S.13 


0.00 


0.00 


0.36 




3.62(R) 



''If the higher M, value of 1.4x10** Mq for NGC 224 l|Bender et alJboOSh and the higher velocity dispersion of 
200 kms~^ for NGC 2787 (see TablefSjl arc used, then the fits are altered as follows: 

23gal(lsqfit): a=3.71 ± 0.37; d=S.15 ± 0.07; ey,™ = 0.29; Ofy = 0.33; a^ax = 3.73 ± 0.38{AI); Omin = 3.66 ± 
0.37(RB) 

8B(lsqfit): a=3.10±0.96; d=8.07±0.17; eyi„ = 0.40; cr/y = 0.45; a^ax = 3.39±0.98{RI); a„i„ = 3.10±0.96{AB) 
23gal(robust): a=3.49; d=8.20; cr/y = 0.33; a^i„ = 3.46(R) 
8B(robust): a=3.66; d=8.28; Ufy = 0.47; a„,ax = 3.74(R) 

Note. — The columns are as follows: (1) galaxy sample, (2) fitting algorithm, (3) slope (a), (4) zero-point (d), 
(5) intrinsic uncertainty in y = logM, {tyin) selected such that Xr = (6) uncertainty in X = logc {ex), 
(7) scatter of points around fit relationship (afy), (S) maximum slope for sample with {weightcent method) and 
(9) minimum slope for sample with (weightcent method). The weightcent method designations are A:Avgerr or 
R:Recent and I=INT, 0=OBS, or B =OBS+INT. The OBS fit for the 8B and 23gal samples are omitted, for 
reasons previously discussed. 
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Table 6. Bootstrap Statistics: M,(cr) 



Parameter 


Isqfit 


Isqfit 


Isqfit 


Isqfit 


robust 


robust 


robust 


robust 


sample 


15E 


23gal 


8B 


Cap8 


15E 


23gal 


8B 


Cap8 




991 


1000 


952 


959 


1000 


1000 


1000 


1000 


a 


Mean 


3.77 


3.81 


4.04 


3.80 


3.73 


3.79 


4.20 


3.69 


Median 


3.84 


3.82 


3.75 


4.03 


3.91 


3.85 


3.66 


3.91 


Adev 


0.27 


0.20 


1.01 


0.54 


0.46 


0.29 


1.73 


1.22 


Sdev 


0.38 


0.27 


1.45 


0.98 


0.65 


0.38 


3.90 


8.12 


Skew 


-0.47 


-0.60 


1.43 


-3.04 


-0.79 


-0.51 


12.74 


-8.30 


Kurt 


4.43 


3.04 


6.41 


13.24 


1.80 


1.89 


247.91 


453.99 


d 


Mean 


8.22 


8.18 


8.18 


8.20 


8.20 


8.18 


8.21 


8.12 


Median 


8.22 


8.18 


8.20 


8.22 


8.16 


8.16 


8.28 


8.26 


Adev 


0.06 


0.05 


0.12 


0.11 


0.09 


0.06 


0.21 


0.24 


Sdev 


0.08 


0.06 


0.15 


0.13 


0.11 


0.08 


0.44 


1.74 


Skew 


-0.54 


-0.17 


-0.40 


-0.79 


0.54 


0.59 


12.27 


-19.22 


Kurt 


1.44 


0.24 


0.13 


0.94 


0.21 


0.60 


224.60 


388.43 




Mean 


0.31 


0.35 


0.40 


0.39 


0.25 


0.28 


0.33 


0.37 


Median 


0.29 


0.35 


0.38 


0.36 


0.24 


0.27 


0.30 


0.30 


Adev 


0.06 


0.05 


0.10 


0.09 


0.04 


0.04 


0.10 


0.15 


Sdev 


0.07 


0.06 


0.12 


0.12 


0.06 


0.05 


0.27 


1.05 


Skew 


0.95 


0.33 


0.53 


1.29 


0.48 


0.09 


16.77 


18.61 


Kurt 


0.76 


-0.02 


0.43 


3.23 


0.11 


-0.03 


339.89 


357.76 




Mean 


0.25 


0.30 


0.32 


0.32 










Median 


0.24 


0.30 


0.30 


0.30 










Adev 


0.07 


0.06 


0.11 


0.12 










Sdev 


0.09 


0.07 


0.14 


0.15 










Skew 


0.65 


0.06 


0.44 


0.77 










Kurt 


0.16 


-0.10 


0.09 


1.24 











Note. — The results from a bootstrap (1000 runs) for each galaxy sample for the M,(cr) 
predictor using the default weightcent Isqfit and robust fitting algorithms arc given. The 
mean, median, average deviation (Adev), standard deviation (Sdev), skew and kurtosis 
(Kurt) are given for the slope (a), zero-point (d), calculated scatter (cfy) and intrinsic un- 
certainty in logM, (eyin) needed to obtain a = 1-0 the least-squares fitting (counting 
only unique galaxies in the DOF). 
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Table 7. Fitting Parameters: log M, = ai log j^ + a2log^ + log ^ + d 



Sample 


ai 


(12 


0-3 


d 






Fay 


Sig 


15E 


O.OOiO.OO 


0.05±0.24 


3.73±0.71 


8.18±0.08 


0.25 


0.29 


0.92 


NSS 


23gal 


O.OOiO.OO 


0.28±0.16 


3.16it0.46 


8.15±0.06 


0.23 


0.30 


1.01 


NSS 


8B 


O.OOibO.OO 


0.94±0.33 


2.90di0.52 


8.32±0.14 


0.16 


0.24 


2.42 


75SS 


15E 


0.03±0.18 


O.OOitO.OO 


3.90it0.50 


8.19it0.09 


0.25 


0.29 


0.93 


NSS 


23gal 


-0.16±0.12 


O.OOitO.OO 


3.56it0.35 


8.15it0.06 


0.24 


0.30 


0.96 


NSS 


8B 


-1.02±0.24 


O.OOitO.OO 


4.68it0.61 


8.54±0.1.5 


0.10 


0.22 


3.06 


75SS 


15E 


-0.71±0.36 


O.OOitO.OO 


O.OOiO.OO 


8.20±0.20 


0.68 


0.70 


0.16 


99SS 


23gal 


-0.62±0.29 


O.OOitO.OO 


O.OOiO.OO 


8.13±0.15 


0.70 


0.72 


0.17 


99SS 


8B 


0.68±0.77 


O.OOitO.OO 


O.OOitO.OO 


7.59it0.33 


0.66 


0.71 


0.28 


90SS 


15E 


O.OOitO.OO 


1.07it0.25 


O.OOitO.OO 


8.13it0.15 


0.49 


0.53 


0.29 


95SS 


23gal 


O.OOitO.OO 


1.13it0.21 


O.OOitO.OO 


S.lOitO.ll 


0.49 


0.53 


0.32 


99SS 


SB 


O.OOitO.OO 


1.56it0.90 


O.OOitO.OO 


8.20it0.31 


0.56 


0.59 


0.41 


75SS 


15E 


2.48±0.88 


3.31itl.l6 


-0.25itl.48 


8.23it0.07 


0.17 


0.24 


1.38 


NSS 


23gal 


0.77±0.51 


1.32±0.70 


1.91±0.93 


8.13it0.06 


0.22 


0.28 


1.17 


NSS 


SB 


-1.33±0.68 


-0.32±0.67 


5.32±1.37 


8.60±0.20 


0.12 


0.24 


2.47 


75SS 


15E 


2.34±0.32 


3.12it0.29 


O.OOitO.OO 


8.23it0.07 


0.16 


0.23 


1.50 


75SS 


23gal 


1.72±0.26 


2.68it0.26 


O.OOitO.OO 


8.11it0.06 


0.24 


0.29 


1.05 


NSS 


8B 


1.47±0.42 


2.39it0.60 


O.OOitO.OO 


7.98it0.20 


0.30 


0.34 


1.21 


NSS 



Note. — Exploration of multivariate fitting parameter combinations. The columns are as follows: 
(1) galaxy sample, (2) ai, (3) a2, (4) as, (5) d, (6) eyi„ selected such that Xr = 1-0, (7) scatter of 
points around fit relationship (afy), (8) F-tcst ratio based on comparison of Af, (a) and M, (fit) and 
(9) the F-tcst probabilistic comparison that the fits are different: NSS (not significantly different), 
75SS (different at 75% level), 90SS (different at 90% level), 95SS (different at 95% level) and 99SS 
(different at 99% level). All fits are determined using the Isqfit algorithm with default weightcent 
selections. 
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Table 8. 


, Fittin} 


5 Parameters: log 


M, = ai lo^ 


' Xo 


+ (32 log 


r Re 
' ReO 


+ d 


X 


Sample 


ai 


a2 


d 






Fay 


Sig 


IlRl 
IlRl 
IlRl 


15E 
23gal 
8B 
Caps 


0.92±0.07 
0.91±0.0S 
0.76±0.16 
0.95±0.13 


O.OOiO.OO 
O.OOiO.OO 
O.OOiO.OO 
O.OOiO.OO 


S.22±0.06 
S.14±0.06 
7.96±0.12 
S.25±0.12 


o.is 

0.24 
0.27 
0.27 


0.25 
0.28 
0.32 
0.32 


1.28 
1.10 
1.39 
1.26 


NSS 
NSS 
NSS 
NSS 




15E 
Cap8 


0.64±0.10 
0.65±0.17 


O.OOiO.OO 
O.OOiO.OO 


S.23±0.10 
S.24±0.19 


0.35 
0.48 


0.39 
0.51 


0.52 
0.51 


75SS 
75SS 


-f 2 r2 n3 

^ best^e^e 
Y2 r2p3 


15E 
23gal 
8B 
Caps 


0.63±0.06 
0.60±0.06 
0.47±0.15 
0.65±0.09 


O.OOiO.OO 
O.OOiO.OO 
O.OOiO.OO 
O.OOiO.OO 


S.1S±0.07 
S.16±0.07 
S.03±0.17 
S.23±0.12 


0.22 
0.29 
0.39 
0.26 


0.29 
0.34 
0.44 
0.32 


0.97 
0.78 
0.73 
1.31 


NSS 
NSS 
NSS 
NSS 


leRe 

leRe 
leRl 


15E 
23gal 
SB 
Caps 


1.07±0.14 
1.13±0.11 
1.24±0.2S 
1.02±0.22 


O.OOiO.OO 
O.OOiO.OO 
O.OOiO.OO 
O.OOiO.OO 


S.17±0.09 
S.13±0.07 
S.12±0.14 
S.15±0.16 


0.30 
0.28 
0.28 
0.41 


0.35 
0.34 
0.33 
0.45 


0.65 
0.79 
1.35 
0.65 


75SS 
NSS 
NSS 
NSS 


'^l JlRl 

best e e 
-v-2 r2 n3 
^best^e^e 
•v-2 1-2 n3 

best-'e-n-e 
•Y-2 r2 n3 


15E 
23gal 
SB 
Caps 


0.73±0.13 
0.53±0.10 
0.41±0.16 
0.79±0.19 


-0.24±0.26 
0.19±0.22 
0.74±0.71 
-0.29±0.34 


S.21±0.0S 
S.15±0.07 
S.20±0.24 
S.27±0.13 


0.22 
0.29 
0.39 
0.27 


0.29 
0.34 
0.43 
0.32 


0.95 
0.76 
0.76 
1.26 


NSS 
NSS 
NSS 
NSS 



Note. — The columns are as follows: (1) fit parameter, (2) galaxy sample, (3) ai, (4) 02, (5) d, 
(6) ej/in, (7) o'fy, (8) the F-test ratio, and its (9) significance as in Table[7] All fits are determined 
using the Isqfit algorithm with default weightcent. 
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Table 9. Fitting Parameters: logM, = a log -et- + d 



Sample 


Algorithm 


a 


d 










Fay 


Sig 


15E 


Isqfit 


0.62±0.06 


8.13±0.07 


0.21 


0.28 


0.62±0.07(AI) 


0.59±0.04{RO) 


1.01 


NSS 


23gaP 


Isqfit 


0.58±0.05 


8.13±0.07 


0.28 


0.32 


0.59±0.06(AI) 


0.57±0.05(RB) 


0.84 


NSS 




Isqfit 


0.47±0.15 


8.05±0.17 


0.38 


0.43 


0.50±0.14(AI) 


0.47±0.15(AB) 


0.78 


NSS 


Caps 


Isqfit 


0.62±0.09 


8.16±0.12 


0.27 


0.32 


0.62±0.09(AB) 


0.61±0.09{RI) 


1.24 


NSS 


15E 


robust 


0.59 


8.14 


0.0 


0.28 




0.59{R) 


1.0 


NSS 


23gal'' 


robust 


0.59 


8.12 


0.0 


0.33 




0.58(R) 


0.83 


NSS 


8B=- 


robust 


0.66 


8.07 


0.0 


0.48 




0.65(R) 


0.59 


NSS 


Cap8 


robust 


0.59 


8.14 


0.0 


0.33 


0.60(R) 




1.20 


NSS 



''If the higher M. value of 1.4xlO*M0 for NGC 224 l|Bender et al.ll2005h and the higher velocity dispersion of 200 
kms~^ for NGC 2787 (see Table|3]l are used, then the fits are altered as follows: 

23gal(lsqfit): a=0.58 ± 0.05; d=8.15 ± 0.07; eyi„ = 0.27; afy = 0.31; amax = 0.58 ± 0.06(AI); a„i„ = 0.57 ± 0.05(RB) 
8B(lsqfit): a=0.49 ± 0.14; d=8.11 ± 0.16; eyi„ = 0.37; afy = 0.40; amax = 0.53 ± 0.13(RI); a™„ = 0.49 ± 0.14(AB) 
23gal(robust): a=0.59; d=8.14; afy = 0.31; a^in = 0.59{R) 
8B(robust): a=0.71; d=8.09; afy = 0.49; a^in = 0.69{R) 

Note. — The columns are as follows: (1) galaxy sample, (2) fitting algorithm, (3) slope (a), (4) zero-point (d), (5) 
Eyini (6) cTfy, (7) maximum slope for sample with (weightcent) , (8) minimum slope for sample with (weightcent) , (9) 
F-test ratio and its (10) significance specified as in Table [7] 



-57- 



Table 10. Fitting Parameters: logM, oc alog 



Sample 




Cx =0.1 






C,r = 0.3 






a 






a 








X=E, 


15E 


0.62(0.62,0.60) 


0.20 


0.28 


0.65(0.66,0.64) 


0.10 


0.28 


23gal 


0.58(0.59,0.58) 


0.27 


0.32 


0.61(0.62,0.61) 


0.21 


0.33 


8B 


0.48(0.50,0.48) 


0.38 


0.43 


0.51(0.53.0.51) 


0.35 


0.43 


Cap8 


0.63(0.6.3,0.62) 


0.26 


0.32 


0.65(0.66,0.64) 


0.19 


0.33 


X=M{10Re) 


15E 


0.83(0.84,0.83) 


0.27 


0.34 


0.92(0.94,0.91) 


0.11 


0.35 


23gal 


0.85(0.86,0.84) 


0.27 


0.34 


0.95(0.97,0.95) 


0.11 


0.35 


8B 


0.90(0.94,0.90) 


0.34 


0.40 


1.14(1.14,1.14) 


0.18 


0.42 


Caps 


0.82(0.83,0.81) 


0.34 


0.40 


0.88(0.91,0.87) 


0.24 


0.40 


X=$(ile/8) 


15E 


■•• (2.23,2.16) 






••• (2.27,2.27) 






23gal 


1.69(1.69,1.66) 


.32 


0.39 


... ( ... ) 






8B 


1.01(1.09,1.01) 


0.42 


0.47 


1.45(1.49,1.45) 


0.25 


0.52 


Cap8 


2.35(2.37,2.20) 


0.11 


0.32 


••• (2.45,2.42) 







Note. — The columns are as follows: (1) galaxy sample; (2) the slope (max, min), (3) 
requisite eyi„ for Xr = 1.0, (4) scatter {Tfy) for ex = 0.1; the (5) slope(max,min), (6) eyi„, 
(7) afy for ex = 0.3. 
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Table 11. Bootstrap Statistics: Eg, $(-Re/8), and M(10i?e) with Isqfit 



Parameter 


Eg 


Eg 


Eg 


Eg 










M 


M 


M 


M 


sample 


15E 


23gal 


8B 


Cap8 


15E 


23gal 


8B 


Cap8 


15E 


23gal 


8B 


Cap8 




998 


1000 


984 


927 


979 


1000 


983 


839 


1000 


1000 


958 


958 


a 


Mean 


0.62 


0.57 


0.42 


0.60 


2.03 


1.55 


0.88 


2.13 


0.83 


0.82 


0.75 


0.76 


Median 


0.61 


0.57 


0.42 


0.63 


2.07 


1.54 


0.90 


2.21 


0.82 


0.82 


0.74 


0.83 


Adev 


0.04 


0.05 


0.12 


0.08 


0.20 


0.20 


0.30 


0.25 


0.07 


0.06 


0.21 


0.16 


Sdev 


0.06 


0.07 


0.15 


0.15 


0.25 


0.25 


0.38 


0.35 


0.10 


0.08 


0.26 


0.31 


Skew 


1.12 


0.00 


0.19 


-2.43 


-0.47 


0.13 


0.25 


-1.77 


0.65 


0.15 


-0.03 


-3.61 


Kurt 


4.45 


1.42 


0.45 


12.14 


-0.03 


0.05 


2.40 


5.21 


6.77 


1.04 


0.19 


15.51 


d 


Mean 


8.15 


8.15 


8.08 


8.17 


8.21 


8.19 


8.01 


8.28 


8.12 


8.10 


8.09 


8.15 


Median 


8.16 


8.15 


8.09 


8.19 


8.20 


8.20 


8.02 


8.29 


8.12 


8.11 


8.10 


8.16 


Adev 


0.06 


0.06 


0.10 


0.09 


0.06 


0.07 


0.11 


0.08 


0.07 


0.05 


0.08 


0.10 


Sdev 


0.08 


0.07 


0.12 


0.11 


0.08 


0.09 


0.14 


0.10 


0.10 


0.07 


0.10 


0.13 


Skew 


-1.03 


-0.01 


-0.34 


-0.88 


0.31 


-0.06 


-0.68 


0.19 


-1.14 


-0.10 


-0.59 


-0.63 


Kurt 


2.88 


0.08 


0.08 


0.84 


-0.45 


-0.03 


2.39 


0.24 


4.43 


-0.12 


1.11 


0.88 




Mean 


0.32 


0.38 


0.47 


0.38 


0.30 


0.49 


0.51 


0.38 


0.39 


0.39 


0.44 


0.42 


Median 


0.32 


0.38 


0.46 


0.37 


0.29 


0.48 


0.51 


0.38 


0.38 


0.39 


0.43 


0.41 


Adev 


0.04 


0.04 


0.10 


0.07 


0.05 


0.07 


0.10 


0.08 


0.05 


0.04 


0.09 


0.10 


Sdev 


0.05 


0.05 


0.13 


0.09 


0.06 


0.08 


0.13 


0.10 


0.07 


0.05 


0.12 


0.13 


Skew 


0.56 


0.29 


0.53 


0.88 


0.46 


0.29 


0.36 


1.03 


0.35 


0.34 


0.58 


0.89 


Kurt 


0.68 


0.19 


1.50 


3.14 


0.51 


0.29 


1.13 


3.16 


0.32 


0.10 


0.50 


1.19 




Mean 


0.26 


0.34 


0.41 


0.30 


0.23 


0.45 


0.45 


0.32 


0.33 


0.34 


0.38 


0.36 


Median 


0.26 


0.33 


0.40 


0.31 


0.23 


0.45 


0.45 


0.32 


0.33 


0.33 


0.37 


0.35 


Adev 


0.05 


0.05 


0.11 


0.09 


0.07 


0.07 


0.11 


0.10 


0.06 


0.05 


0.11 


0.13 


Sdev 


0.06 


0.06 


0.15 


0.12 


0.08 


0.09 


0.15 


0.13 


0.08 


0.06 


0.14 


0.16 


Skew 


0.23 


0.17 


0.41 


0.41 


-0.03 


0.10 


0.21 


0.18 


0.18 


0.24 


0.35 


0.51 


Kurt 


0.38 


0.14 


0.94 


1.31 


0.00 


0.28 


0.77 


1.26 


0.14 


0.10 


0.20 


0.29 



Note. — For Isgfit (default weightcent) the results are given from 1000 bootstrap 'runs' for e3x:h galaxy sample 
when fitting SMBH mass against the gravitational binding energy, gravitational potential at Re/8 and the bulge 

mass, Af(lO-Re). The mean, median, average deviation (Adev), standard deviation (Sdev), skew and kurtosis 
(Kurt) are given for the slope (a), zero-point (d), scatter {cfy) and intrinsic uncertainty (eyin) needed to obtain 
a Xr = 1-0 foi' the least-squares fitting (counting only unique galaxies in the DOF). 
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Table 12. Bootstrap Statistics: Eg, $(i?e/8), and M(10i?e) with robust 



Parameter 




Eg 


Eg 


E, 


$ 


4> 






M 


M 


M 


M 


sample 


15E 


23gal 


8B 


Caps 


15E 


23gal 


8B 


Cap8 


15E 


23gal 


SB 


Cap8 


a 


Mean 


0.60 


0.57 


0.43 


0.68 


2.03 


1.63 


0.79 


2.15 


0.81 


0.80 


0.82 


0.80 


Median 


0.59 


0.58 


0.42 


0.62 


2.13 


1.76 


0.94 


2.48 


0.77 


0.77 


0.77 


0.89 


Adev 


0.07 


0.08 


0.25 


0.24 


0.30 


0.33 


0.67 


0.55 


0.12 


0.12 


0.46 


0.17 


Sdev 


0.10 


0.11 


0.35 


3.06 


0.37 


0.39 


2.59 


1.85 


0.29 


0.16 


1.43 


0.38 


Skew 


0.85 


-0.14 


4.26 


29.85 


-0.48 


-0.55 


-17.84 


-17.92 


15.56 


0.27 


8.72 


1.67 


Kurt 


5.16 


1.31 


54.38 


928.65 


0.64 


-0.64 


353.75 


538.13 


380.59 


0.33 


373.24 


59.03 


d 


Mean 


8.16 


8.14 


8.09 


8.18 


8.18 


8.18 


7.94 


8.29 


8.11 


8.11 


8.10 


8.18 


Median 


8.14 


8.12 


8.07 


8.21 


8.20 


8.15 


8.03 


8.27 


8.05 


8.07 


8.07 


8.20 


Adev 


0.10 


0.08 


0.20 


0.23 


0.09 


0.09 


0.30 


0.12 


0.14 


0.11 


0.24 


0.15 


Sdev 


0.14 


0.11 


0.45 


1.46 


0.11 


0.11 


1.77 


0.86 


0.27 


0.12 


1.18 


0.18 


Skew 


-0.95 


-0.16 


13.75 


-9.20 


0.28 


0.50 


-19.13 


24.36 


-14.74 


0.19 


0.98 


-0.86 


Ktu-t 


2.72 


2.11 


260.07 


332.33 


0.51 


-0.63 


.388.99 


657.81 


351.03 


-0.22 


362.95 


12.06 




Mean 


0.26 


0.31 


0.40 


0.34 


0.24 


0.39 


0.46 


0.32 


0.32 


0.32 


0.40 


0.32 


Median 


0.26 


0.31 


0.40 


0.29 


0.24 


0.39 


0.43 


0.30 


0.32 


0.32 


0.38 


0.32 


Adev 


0.04 


0.03 


0.11 


0.13 


0.05 


0.05 


0.15 


0.11 


0.05 


0.03 


0.15 


0.09 


Sdev 


0.05 


0.04 


0.20 


0.90 


0.06 


0.07 


0.86 


0.53 


0.11 


0.04 


0.76 


0.11 


Skew 


0.49 


0.69 


14.34 


17.83 


0.17 


0.03 


22.45 


23.68 


18.86 


0.31 


20.08 


0.37 


Kurt 


0.58 


3.34 


346.76 


320.34 


-0.18 


-0.08 


547.67 


632.84 


496.01 


0.25 


424.37 


0.30 



Note. — For robust the results are given from 1000 bootstrap 'runs' for each galaxy sample when fitting SMBH mass 
against the gravitational binding energy, gravitational potential at i?e/8 and the bulge mass, M(lOile). The mean, median, 
average deviation (Adev), standard deviation (Sdev), skew and kurtosis (Kurt) are given for the slope (a), zero-point (d), 
and scatter {<Tfy)- 
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Table 13. 


Fitting 


Parameters: log 


;M. = 


aloe 


+ d 








Sample 


Algorithm 


a 


d 










Fay 


Sig 


$(20pc) 
3>(20pc) 
$(20pc) 
<I>(20pc) 


15E 
23gal 
8B 
Cap8 


Isqfit 
Isqfit 
Isqfit 
Isqfit 


1.85±0.16 
1.49±0.17 
0.97±0.32 
1.84±0.27 


8.22±0.07 
8.22±0.08 
8.02±0.18 
8.23±0.13 


0.19 
0.33 
0.41 
0.29 


0.26 
0.37 
0.45 
0.34 


1.88±0.18(AI) 
1.50±0.17{AI) 
1.04±0.31(RI) 
1.84±0.27{AB) 


1.74±0.13(RO) 
1.48±0.17{RB) 
0.97±0.32(AB) 
1.80±0.19(RO) 


1.22 
0.66 
0.70 
1.14 


NSS 
75SS 
NSS 
NSS 


<I>(100pc) 
*(100pc) 
*(100pc) 
$(100pc) 


15E 
23gal 
8B 
Cap8 


Isqfit 
Isqfit 
Isqfit 
Isqfit 


1.48±0.14 
1.34±0.14 
0.96±0.32 
1.44±0.24 


8.24±0.07 
8.21±0.08 
8.02±0.18 
8.24±0.14 


0.22 
0.31 
0.41 
0.32 


0.28 
0.35 
0.45 
0.36 


1.51±0.16{AI) 
1.36±0.15(AI) 
1.03±0.31(RI) 
1.44±0.24(AB) 


1.38±0.11(RO) 
1.33±0.14(RB) 
0.96±0.32(AB) 
1.41±0.25(RI) 


1.02 
0.72 
0.71 
0.99 


NSS 
75SS 
NSS 
NSS 


HRe) 


15E 
23gal 
8B 
Cap8 


Isqfit 
Isqfit 
Isqfit 
Isqfit 


1.95±0.20 
1.53±0.21 
0.94±0.37 
2.03±0.34 


8.26±0.08 
8.19±0.09 
7.92±0.18 
8.31±0.14 


0.24 
0.40 
0.46 
0.32 


0.29 
0.43 
0.50 
0.37 


2.09±0.17(AO) 
1.53±0.21(AB) 
1.00±0.37(RI) 
2.30±0.30(AO) 


1.84±0.23(RI) 
1.50±0.22(RI) 
0.94±0.37(AB) 
1.87±0.41(RI) 


0.92 
0.49 
0.56 
0.97 


NSS 
90SS 
75SS 
NSS 


-I-CRe/S) 
^{Re/S) 
*(fie/8) 


15E 
23gal 
8B 
Cap8 


Isqfit 
Isqfit 
Isqfit 
Isqfit 


2.09±0.17 
1.60±0.20 
0.98±0.35 
2.22±0.30 


8.19±0.06 
8.18±0.08 
7.96±0.18 
8.26±0.12 


0.18 
0.36 
0.43 
0.25 


0.25 
0.39 
0.48 
0.32 


2.14±0.14(AO) 
1.60±0.20{AI) 
1.05±0.35(RI) 
2.37±0.23(AO) 


2.00±0.21(RI) 
1.57±0.20{RB) 
0.98±0.35(AB) 
2.06±0.38(RI) 


1.29 
0.58 
0.63 
1.31 


NSS 
75SS 
NSS 
NSS 


*(fie/8) 
^(iJe/S) 
^(Re/S) 
*(Re/8) 


15E 
23gal 
8B 
Cap8 


robust 
robust 
robust 
robust 


2.14 
1.79 
1.37 
2.47 


8.17 
8.10 
8.00 
8.26 


0.00 
0.00 
0.00 
0.00 


0.25 
0.41 
0.51 
0.34 


2.13(R) 
1.84(R) 


1.36(R) 
2.47(R) 


1.25 
0.54 
0.53 
1.12 


NSS 
90SS 
75SS 
NSS 



Note. — 



Gravitational Potential at radii 'X', as specified. See Table |9]for a description of columns 2-11. 
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Table 14. 


Fitting '. 


Parameters: 


: log 


;M. = 








X 


Sample 


Algorithm 


a 


d 










Fay 


Sig 


M(R^) 

M{Re) 

M{R,) 


15E 
23gal 
8B 
Caps 


Isqfit 
Isqfit 
Isqfit 
Isqfit 


O.S2±0.11 
O.SO±0.12 
O.S3±0.25 
0.79±0.16 


S.06±0.10 
S.12±0.09 
S.09±0.1S 
S.05±0.16 


0.31 
0.41 
0.3S 
0.39 


0.36 
0.44 
0.42 
0.43 


O.S3±0.10(RI) 
0.S1±0.12(AI) 
O.SS±0.24(AI) 
O.SO±0.10(RO) 


0.77±0.0S(RO) 
O.SO±0.11(RB) 
O.S3±0.26{RB) 
0.79±0.17{AI) 


0.63 
0.45 
O.Sl 
0.69 


75SS 
95SS 
NSS 
NSS 


M(10i?e) 
M{10Re) 

M{10Re) 
M(10R,) 


15E 
23gal 
SB 
Caps 


Isqfit 
Isqfit 
Isqfit 
Isqfit 


O.S2±0.10 
0.S3±0.0S 
O.SS±0.25 
0.S1±0.15 


S.09±0.09 
S.0S±0.07 
S.09±0.17 
S.09±0.14 


0.2S 
0.29 
0.35 
0.35 


0.34 
0.34 
0.40 
0.40 


O.S2±0.10{RI) 
O.S5±0.09(AI) 
0.93±0.23(AI) 
O.S1±0.09{RO) 


0.79±0.07(RO) 
0.S3±0.0S{RB) 
O.SS±0.25{AB) 
O.SO±0.15(RI) 


0.70 
0.76 
0.91 
0.S3 


NSS 
NSS 
NSS 
NSS 


M(10Re) 
MiWRe) 
MilORe) 
M{10Re) 


15E 
23gal 
SB 
Caps 


robust 
robust 
robust 
robust 


0.77 
0.77 
1.23 
0.77 


S.04 
S.05 
S.07 
S.05 


0.00 
0.00 
0.00 
0.00 


0.35 
0.35 
0.47 
0.40 


0.7S(R) 
0.7S{R) 


0.76{R) 
1.19(R) 


0.64 
0.73 
0.62 
O.Sl 


75SS 
75SS 
NSS 
NSS 


Af(lOfcpc) 
M(lOfcpc) 
M(lOfcpc) 
M(lOfepc) 


15E 
23gal 
SB 
Caps 


Isqfit 
Isqfit 
Isqfit 
Isqfit 


0.92±0.11 
0.94±0.10 
0.91±0.26 
0.90±0.17 


S.15±0.09 
S.12±0.07 
S.05±0.16 
S.15±0.15 


0.29 
0.30 
0.36 
0.36 


0.34 
0.35 
0.40 
0.41 


0.93±0.12(AI) 
0.96±0.10(AI) 
0.96±0.25(AI) 
0.90±0.17(AB) 


0.S5±0.0S(RO) 
0.93±0.09{RB) 
0.91±0.26{AB) 
O.S9±0.17(RI) 


0.69 
0.74 
O.SS 
0.7S 


NSS 
NSS 
NSS 
NSS 


M{kpc) 
M(kpc) 
M(kpc) 
M(kpc) 


15E 

23gal 

SB 

Caps 


Isqfit 
Isqfit 
Isqfit 
Isqfit 


1.35±0.15 
1.29±0.15 
0.97±0.33 
1.26±0.25 


S.24±0.0S 
S.16±0.0S 
7.94±0.17 
S.20±0.16 


0.27 
0.34 
0.41 
0.3S 


0.32 
0.3S 
0.46 
0.42 


1.37±0.17(AI) 
1.31±0.16(AI) 
1.04±0.32(RI) 
1.26±0.25(AB) 


1.23±0.11(RO) 
1.27±0.15{RB) 
0.97±0.33{AB) 
1.23±0.27(RI) 


0.77 
0.62 
0.69 
0.74 


NSS 
75SS 
NSS 
NSS 



Note. — 



Mass enclosed by spherical radius X, as specified in column 1. See Table[9]for a description of columns 2-11. 
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Table 15. Maximum and Minimum Limits for Galaxy Parameters: Inclination 



Galaxy 


e 


Rcf 


Smax 


Ref 


Smin 


Ref 


MAX/MIN e 


NGC 221 


70 


1 






50 


2 


■■■ /28.6% 


NGC 224 


77 


3,4,5 


78 


6 


55 




1.3% /28.6% 


NGC 1023 


90 


7 






60.4 


8 


■■• /32.9% 


NGC 1068 


38.3 


8 


46 


9 


25 


10 


20.1% /34.7% 


NGC 3379 


90 


2,11 






60*= 


12 


■■• /33.3% 


NGC 3384 


90 


13 






64 


14 


■■• /28.9% 


NGC 4261 


64 


15 






30 


12 


■■• /53.1% 


NGC 4342 


90 


16 






78^^ 


17 


■■• /13.3% 


NGC 4459 


41 


18 


47 


2,18b 






14.6% / ■ ■ ■ 


NGC 4486 


42 


19 


SI'* 


20 






21.4% / ■ ■ ■ 


NGC 6251 


76 


21 






68 


22 


■■• /10.5% 


NGC 7457 


90 


13 






64 


2 


■■■ /28.9% 


IC 1459 


90 


23 






60 


24 


■■• /33.3% 



'^Nuclear Disk 

^Capaccioli et aT 1 1 I1991I) gives a value of 31.5°based on estimate of disk. 
' IScorza &: van den Bosc 3 lll998l ') gives 83° for the nuclear disk. 
' ^Cappellari et ah I ll2006h gives a value of 90° . 



Note. — Description of columns: (1) galaxy name; (2) incli- 
nation (in °), from (3); (4) maximum inclination (in °), from (5); 
(6) minimum inclination (in °), from (7); and (8) percentage variation 

[ioo{9max - eoBs)/{eoBs)]/[ioo{eoBs - eMiN)/(eoBs)]- 



References. — (1) Iverolme et ahl ll2002h (2) ICappellari et ahl l l20o3) . (3) 
KennicuttI lll987l'). (4) iTremainel fl995f ) . (51 iBacon et ahl (2 001), (6) 
(7) Bower et ah i 200il). (8*) ISanchez-Portal et ah. (.2000) 



Walterbos^^^_^^^^ 

Marcum et ahl i 200lll 




T^kamiva fc Sofuj 1I2OOOI). (10*) ICreenhill et ahl lll996l) , (11*) Gebhardt et al 



(1 2) Ivan der Marel Binnev fc Daviej h990h. (13) 
(14) iBusarello et al.Nl996h . (15 



Ferrarese Jord, h Jaffd ||1996|) 



ScorzeTfc'rair'den^oschr (1991 



Gebhardt et al, 
1 



Cretton & van den Bosch' (19 991), ( 1 7) . 

^arzi etah (2001), fl8b ) Sarziet ahl ll200ll) mo deled gas disk,7l 9) Harms et ahl 
Tl994h. (20) iMacchetto et ah! lll997^. (21) iFerrarese fc Ford] (Il999l). (22) 
ICrane fc Vernej lll997h . (23) ICappellari et ahl (I2OO2I) . (24) IVerdoes Kleiin et ahl 
l l2000h . 
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Table 16. Maximum and Minimum Limits for Galaxy Parameters: Axis Ratio 



Galaxy q Ref Qmax Ref qMiN Ref MAX/MIN q 



NGC 


221 


0.73 


1 


0.87 


NGC 


1023 


0.72 


4 


0.799 


NGC 


2778 


0.78 


7 


0.889 


NGC 


3115 


0.60 


9 




NGC 


3377 


0.50 


7 


0.64 


NGC 


3379 


0.90 


10 


1.0 


NGC 


3384 


0.60 


7 


0.82 


NGC 


4261 


0.75 


11 


0.89 


NGC 


4291 


0.76 


7 


0.85 


NGC 


4473 


0.62 


7 




NGC 


4564 


0.70 


7 




NGC 


4486 


0.915 


14 


1.0 


NGC 


6251 


0.85 


16 


0.96 


NGC 


7052 


0.55 


2 


0.89 


NGC 


7457 


0.70 


7 


0.936 



2 


0.63 


3 


19.2% /13.7% 


5 


0.649 


6 


11.0% /9.9% 


6 


0.72 


8 


14.0% /7.7% 




0.452 


5 


■•• /24.7% 


6 


0.491 


5 


28.0% /1. 8% 


3 


0.85 


2 


11.1% /5.6% 


5 


0.585 


6 


36.7% /2.5% 


12 


0.74 


13 


18.7% /1. 3% 


8 


0.74 


6 


11.8% /2.6% 




0.531 


6 


■■■ /14.4% 




0.42 


12 


■•• /40.0% 


3,13 


0.65 


15 


9.3% /29.0% 


17 


0.83 


12 


12.9% /2.4% 


18 


0.44 


18 


61.8% /20.0% 


6 


0.629 


5 


33.7% /10.1% 



Note. — Description of columns: (1) galaxy name; (2) observed 

axis ratio (q), from (3); (4) maximum axis ratio, from (5); (6) 

minimum axis ratio, from (7); and (8) percentage variation given as 
[lOOiQMAx - qoBs)/{qoBs)]/['^oo{qoBs - qMIN)/{qOBs)]- 



References 



1 



I Verolme et al] tOO'j ). 



Iv 

value, (3) iFaber et all lll989l ') ESQ value 
iLauer et al] ll2005h 
e llipticity 



average 
hardt, (8) |r 
(10 
Il98l 



(2) iFaber et al.l ll 19891) FW 

(4) iBower et alT (l200lh . [5) 

average ellipticity r > rj, , (6) iLaueret al.l 1120051) 
< r;,, (7) priv ate communication with K. Geb- 
acon. Monnet. fc Simien l ll985 
1989 



central value of 



Bacon, Monnet, 



Faber et al.. . . ^ ^ 

bimion 119850 de Vaucouleurs 



( 9) iKormendv fc Richstonj lll992l) 
Bender. Dobereiner. fc MoUenhofj 



Capaccioli" fl97! 



van der Marell lll991^. fll) IFerrarese et al.l 1 


19941): IFerrarese, Ford, fc jaffe 


11996 


), (12) Ide Vaucouleurs et al.l il99lh, (13) 


Binggeli, Sandage, fc Tammann 


(1985 


). (14) central value of iFaber et al.l 


(Il989l): Peletier et al.l (Il990l): 


Bacon. Monnet. fc SimienI lll985l), (15) Caon. Capaccioli, fc Rampazzol lll99Cll) SB 



(18) SB profile of lde Juan. Colina. fc Perez-Fournonni994h 



Table 17. Maximum and Minimum Limits for Galaxy Parameters: Mass-to-Light Ratio 



Galaxy 


^bpC 


Ref 




Rcf 


Color 


Rcf 


T 


Ref 


Color 


Ref 


Band 


MAX/MIN T 




i.D 


1 
i 


9 n 
z.u 


z 


u.yoi^Dvj, 1. ivSoi^vi ) 


Ol^DVJ,4^^VlJ 


1 A 
1.4 





l.lOO^Vl } 


A 
4 


T 
1 


ZO.U/0 /iZ.O/0 


NFl^r^ 99/1 


O.U 




0.0 


7 
( 






A 1 
4.1 


Q 







V 


1 /IS nCfl 
iO/0 /io.U/0 


T^C^C* 091 


7 fi 
1.0 





iz. 


z 


u.yy i^Dv ) 


Q 


4.0 


lU 


n fll f^r-f\ 

U.Dii^vr ) 


1 1 
1 1 


V 


£\A ^oy /Of? 00? 
D4.O/0 /oD.0/0 




O.U 


1 9 


0.0 


lU 














V 


ID.U/o / ■ ■ ■ 


INLjV^ lUDo 


1 9 


ioa 


o.V 


lv5 






n 9 
U.z 


1 n 
lU 


U. ( oi^vr ) 


1 A 
14 


V 


iOU/o /oo.o/o 


IN Z i (0 


O.U 




y 










a A 

D.4 









V 


/on n'V^ 


IN V_i V_/ Z J / 


9 


10 


10.4 


ID 














l> 


001.0 /of 


INLj*.^ olio 


D.y 





7 n 


1 n 
lU 






Q 7 
Z. 1 


1 7 
1 / 






V 


1.4/0 / DU.y/o 


TsJPP QQ77 
i\ 00 t i 


9 Q 


Q 

y 


Q A 


9 
z 


U.ODl^DV J 


Q 


9 

Z.O 


lU 






V 


ZZ4.1/0 /lO.o/o 


ivrr^o QQ7Q 


4.0 


lo 


7 

( .0 


ly 


u.oy i^vr J 


9n 


4.0 


91 
Zl 


U.oy i^vr ) 


9n 

ZU 


V 


c;fi /I n 
OD.0/0 /iu.4/0 


ivr^r^ QQQ/1 


9 

Z.O 



y 


9 fl 
Z.O 


D,1U 














V 


1 9 ncfl / 
iz.u/o / ■ ■ ■ 


iNLrO oDUo 


Q T 
0. ( 



y 


D.y 




z 


U.y4l^Dv J 


q 










V 


Q£? CO/ / 

00.0/0 / ■ ■ ■ 


NGC 4261 


4.8 


22 


14.3 


2 


0.99(bv) 


3 










V 


197.9% / ■ ■ ■ 


NGC 4291 


5.5 


9 


12.7 


2 


0.96(bv) 


3 


4.4 


6 






V 


130.9% /20.0% 


NGC 4459 


6.3 


23 


9.3 


16 






3.6 


24 






B 


47.6% /42.9% 


NGC 4473 


6.0 


9 


7.0 


2 


0.96(bv) 


3 


5.2 


10 






V 


16.7% / 13.3% 


NGC 4486 


4.0 


6 


11.0 


23 


0.96(bv) 


3 










V 


175.0% / • ■ ■ 


NGC 4564 


2.0 


9 


10.8 


2 


0.93(bv),1.161(vi) 


3(bv),4(vi) 


1.9 


6 






I 


440.0% /5.0% 


NGC 4596 


1.2 


25 


8.9 


16 














B 


641.7% / • ■ ■ 


NGC 4697 


4.7 


9,10 


7.6 


2 


0.91(bv) 


3 










V 


61.7% / ■ ■ ■ 


NGC 4742 


2.7 


2 










1.1 


10 


1.28(br) 


26 


B 


■■■ /59.3% 


NGC 5845 


5.5 


9,10 


6.5 


2 


1.02(bv) 


3 


4.8 


6 






V 


18.2% /12.7% 


NGC 6251 


9.7 


27 










7.3 


10 


0.61(vr) 


11 


V 


■■• /24.7% 


NGC 7052 


6.3 


6,28 










2.5 


10 


0.70(ri) 


11 


I 


•■■ /60.3% 



Note. — Description of columns: (1) galaxy name; (2) T^pc in Mq/Lq, from (3) reference; (4) Tmax in ^q/Lq, from (5) reference, with (6) 
color correction, from (7) reference; (8) Tjv/jAr in Mq/Lq, from (9) reference, with (10) color correction, from (11) re ference; all in (12) wavelength- 
band; with (13) percentage variation [100(Tmax - Toss)/(Tf"oBs)]/[100(Toss - '^min)/C^OBs)]- ICappellarT e t al. (2006) T-values (Tcap) 
are excluded from consideration here, although, when available, they are superior to those values given in column 2 (see Sections 12. 21 and 14.21 1. For 
those galaxies which are excluded, the previously discussed value (Section I2.2|l was accepted as the correct value. All values have been rescaled 
to a common distance and wavelength-band. The color corrections are coded such th at color X-Y will be given as 'xy'. In the cases where it 
was necessary to convert T to the default waveband, the solar colors were taken from iLivinestonI ||2000|) such that bvQ = 0.650, ubQ = 0.195, 



uvq = 0.845, vtq = 0.54, viq = 0.88, vkQ = 1.49, and the color terms for the individual galaxies were taken from the literature. 
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value of 1.19., (26) IPoulain fc NietQ 119941') . (27) IFerrarese fc FordI lll999l ). (28) Ivan der Marel fc van den BoschI lll998l) . 



Table 18. Fitting Parameters: log M, = ai[log + log ^ + d 



Sample Fx sigx 



X=t7 

15E 0.33 ± 1.62 3.89 ± 0.43 8.17 ± 0.11 0.25 0.29 0.86 ± 1.77(RI) -2.16 ± 1.27(RO) 0.92 NSS 0.92 NSS 

23gal 0.97 ± 1.22 3.90 ± 0.37 8.13 ± 0.08 0.25 0.30 1.00 ± 1.26(AI) 0.78 ± 1.24(RB) 0.99 NSS 0.99 NSS 

8B 1.68 ± 5.39 4.04 ± 2.51 8.11 ± 0.17 0.33 0.40 3.89 ± 5.26{AI) 1.33 ± 5.24(RB) 0.92 NSS 0.92 NSS 

Cap8 1.07 ill 2.75 4.12 ± 0.89 8.12 ± 0.17 0.35 0.39 1.13 ili 2.77(AI) 0.43 ili 2.64(RO) 0.86 NSS 0.86 NSS 

15E 0.04 ± 0.04 0.65 ± 0.07 8.05 ± 0.10 0.21 0.29 0.06 ± 0.03(AO) 0.03 ± 0.04(AI) 0.96 NSS 0.93 NSS 

23gal 0.05 ± 0.04 0.61 ± 0.06 8.06 ± 0.09 0.27 0.32 0.05 ± 0.04{RI) 0.04 ± 0.04{RB) 0.86 NSS 1.0 NSS 

8B 0.08 ± 0.16 0.57 ± 0.24 8.00 ± 0.22 0.42 0.45 0.09 ± 0.15{AI) 0.07 ± 0.16(RB) 0.71 NSS 0.91 NSS 

Cap8 0.04 ± 0.05 0.69 ± 0.12 8.10 ± 0.14 0.28 0.34 0.04 ± 0.05{AO) 0.04 i 0.05(RO) 1.13 NSS 0.89 NSS 

15E -0.21 ± 0.40 2.06 it 0.18 8.22 ± 0.09 0.19 0.26 0.01 ± 0.50(RI) -0.21 ± 0.40(AB) 1.16 NSS 0.92 NSS 

23gal 0.43 ± 0.48 1.71 ± 0.24 8.10 ± 0.12 0.36 0.39 0.43 ± 0.48(AB) 0.39 ± 0.49{RB) 0.58 75SS 1.0 NSS 

8B 0.25 ± 1.06 1.09 ± 0.59 7.91 ih 0.27 0.48 0.51 0.40 ± 0.99(AI) 0.19 ± 1.07(RB) 0.55 NSS 0.89 NSS 

Cap8 -0.10 ifc 0.80 2.18 ili 0.39 8.27 ifc 0.16 0.29 0.35 0.03 ifc 0.98(RI) -0.48 ifc 0.79(RO) 1.09 NSS 0.84 NSS 

X=M(10/?e) 

ISE 0.08 it 0.09 0.88 it 0.11 8.00 it 0.13 028 035 0.15 it 0.07(AO) 0.07 it 0.09(AI) 066 75SS 094 NSS 

23gal 0.08 it 0.08 0.87 it 0.09 8.02 it 0.09 0.29 0.34 0.08 ± 0.08{RI) 0.07 ± 0.08(AI) 0.76 NSS 1.0 NSS 

8B 0.18 it 0.44 1.00 it 0.38 8.05 ± 0.20 0.38 0.42 0.19 it 0.41(AI) 0.15 it 0.45(RB) 0.80 NSS 0.91 NSS 

Cap8 0.08 it 0.13 0.89 it 0.21 8.03 it 0.18 0.37 0.42 0.09 it O.ll(AO) 0.07 it 0.13(AI) 0.73 NSS 0.91 NSS 



Note. — The log-quadratic fits with the columns as follows: (1) galaxy sample, (2) ai, (3) 02, (4) d, (5) eyin, (6) o-fy, (7) maximum ai (weightcent) 



and (8) minimum ai (weightcent), (9) F-tcst ratio and (10) significance from comparison with the log-linear velocity dispersion predictor and (11) 
F-test ratio and (12) significance from comparison with the log-linear fit for X. All fits determined using the Isqfit algorithm with default weightcent. 
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Table 19. Fitting Parameters: log M, = oi log + a2 log + d 



Sample 


ai 


12 


d 






Fay 


siga 


Fi 


sigi 


15E 


0.88±0.12 


-0.62±0.26 


8.16±0.06 


0.16 


0.24 


1.36 


NSS 


1.36 


NSS 


23gal 


0.57±0.10 


0.02±0.24 


8.13±0.07 


0.28 


0.33 


0.80 


NSS 


0.94 


NSS 


8B 


0.41±0.17 


0.59±0.76 


8.18±0.24 


0.40 


0.44 


0.74 


NSS 


0.96 


NSS 


Cap8 


0.92±0.17 


-0.62±0.33 


8.24±0.11 


0.20 


0.27 


1.73 


NSS 


1.40 


NSS 



Note. — The multivariate (Eg, Re) predictive fits with the columns as follows: (1) galaxy sample, 

(2) ai, (3) 02, (4) d, (5) f-yin, (6) cr/y, (7) F-tcst ratio and (8) significance from comparison with the 
log-linear velocity dispersion predictor and (9) F-test ratio and (10) significance from comparison 
with the log-linear fit (for Eg). All fits are determined using the Isqfit algorithm with default 
weightcent. 
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Table 20. 


Fitting Parameters: 


logM, 


= ai lo; 


^2L 
^ Xo 


+ (32 loi 




+ d 




X 


Sample 


ai 


a2 


d 






Fay 


siga 


Fi 


sigi 


HR^) 


15E 
Cap8 


1.55±0.23 
1.62±0.34 


0.38±0.15 
0.42±0.19 


8.19±0.07 
8.25±0.12 


0.18 
0.23 


0.26 
0.30 


1.18 
1.45 


NSS 
NSS 


1.24 
1.52 


NSS 
NSS 


HRe/8) 
*(Re/8) 
*(Re/8) 
*(i?e/8) 


15E 
23gal 
8B 
Cap8 


1.81±0.26 
1.15±0.22 
0.82±0.35 
1.90±0.41 


0.22±0.16 
0.55±0.17 
0.93±0.72 
0.25±0.22 


8.15±0.07 
8.14±0.07 
8.18±0.24 
8.22±0.12 


0.17 
0.29 
0.41 
0.24 


0.25 
0.34 
0.44 
0.31 


1.33 
0.79 
0.73 
1.34 


NSS 
NSS 
NSS 
NSS 


1.0 
1.32 
1.19 
1.07 


NSS 
NSS 
NSS 
NSS 



Note. — The multivariate ($(ni?,e), nf?e) fit with the (1) gravitational potential at radius R and columns 
(2) - (11) as in Table [T51 columns (1) - (10). 
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Table 21. Fitting Parameters: log M. = ai log '"//"^Y + «2 log + d 



Sample ai a2 d eyi„ Ufy F^y sig^ Fi sigi 

15E 1.54±0.23 -1.17±0.35 8.11±0.07 0.18 0.26 1.19 NSS 1.71 75SS 

23gal 1.09±0.19 -0.45±0.30 8.08±0.07 0.28 0.33 0.82 NSS 1.06 NSS 

SB 0.82±0.33 0.26±0.80 8.14±0.23 0.39 0.43 0.78 NSS 0.87 NSS 

Cap8 1.61±0.29 -1.20±0.42 8.17±0.10 0.20 0.27 1.83 NSS 2.19 75SS 



Note. — 



The multivariate (M(10i?e), lORe) fit with columns (1) - (10) as in Table [T9l 
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Fig. 1. — As a function of radius, from left to right, the surface brightness, the deprojected 
luminosity density, and cilogi//cilogr in the V-band along the major axis for the elliptical 
galaxy NGC 6251. 
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Ellipticals: Mbh-E^ Spiral/Lent Bulges; Mbh-E^ All Galaxies: Mbh-E^ 




2 2.5 2 2.5 2 2.5 



log a log a log a 

Fig. 2. — Top three panels illustrate the relationship between the gravitational binding 
energy and the black hole mass, shown with observational uncertainties on the black hole 
mass, and bottom three panels illustrate the relationship between velocity dispersion {a) and 
the black hole mass for the same galaxies. From left to right, the panels show ellipticals only, 
spiral/lenticular bulges only and, finally, all galaxies. Stars indicate elliptical galaxies, with 
filled stars denoting those galaxies in the 15E sample used for the fitting, and the circles 
indicate bulges, with filled circles denoting those galaxies in the 8B sample used for fitting. 
The fit for the combined galaxy sample {23gal) is denoted by a sohd hue, the fit for eUipticals 
only {15E) is denoted by a dotted line, and the fit for the bulges only {8B) is denoted by a 
dashed line, in all panels. 
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Fig. 3. — Illustration of the residuals from the dispersion fit {M,obs — M,pred)', the top panel 
is the 15E elliptical fit, the middle panel is the 8B spiral/lenticular bulge fit and the bottom 
panel is the 23gal sample fit. The left column plots the residual as a function of velocity 
dispersion, the middle column plots the residual as a function of bulge effective radius, and 
the right column plots the residual as a function of V-band /g. The symbols are as in Figure[2l 
The galaxies with \residuals\ > 0.6 in the top panel are NGC 821 {15E,Cap8; filled star) 
and NGC 2778 (open star), in the middle panel are NGC 2778 and NGC 3115 {8B, filled 
circle) and in the bottom panel are NGC 821, NGC 2778 and NGC 3115. 



6739 10 6789 10 6799 10 6789 10 6789 10 
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M-tj; 8B fit M-o-R,-!,; 8B fit l^-R,-!,; SB fit M-o-R,; 8B fit M-cj-l,; 8B fit 
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Fig. 4. — Observed black hole mass plotted against the predicted black hole mass for multi- 
variate fits. The top row is for the fit based on ellipticals {15E), the middle row is for the 
fit based on bulges {8B) and the bottom row is for the fit based on "all" galaxies {23gal). 
From left to right the columns illustrate the M,{a) fit, the M,{Ie, Re,cr) fit, the M,{Ie,Re) 
fit, the M,{Re, a) fit and the M,{Ie, a) fit. The symbols are as in Figure O 




Fig. 5. — Illustration of the relationship between /f-Rg and black hole mass. The left panel 
shows the ellipticals, the central panel shows spiral/lenticular bulges, and the right panel 
shows all galaxies. Lines and symbols are as in Figure O 



Cops 



15E 




Fig. 6. — Top panel illustrates the I^Rl fit, the middle panel the Tl^^I^Rl fit using the Tbpc 
value s from Table El and the bottom panel the Tl^^^I^R^ fit using the revised Tcap values 
from ICappellari et al.l (120061 ). where available. The left column is for the Cap8 elliptical 
galaxy (sub)sample and the right column is for the 15E elliptical galaxy sample. The lines 
illustrate the best-fit relationship for the galaxies in each panel. The symbols are as in 
Figure [2l 
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Ellipticals: Mbln-T^l/R^.^ 
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Fig. 7. — Same as figure [5], but for T'^I^Rl 
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Ellipticals: Mbh-lj.R/ 
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Fig. 8 



— Same as figure [5], but for leRl 
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Fig. 9. — Residuals from the gravitational binding energy fits {M,obs — M,pred)- The top row 
shows the 15E elliptical fit, the middle row the 8B bulge fit and the bottom row the 23gal 
fit. The left column plots this residual as a function of Eg and the right column plots the 
residual as a function of bulge effective radius. The symbols are as in Figure [2l The galaxies 
with \residuals\ > 0.6 in the top panel {15E fit) are NGC 1068 (open circle), NGC 2778 
(open star), NGC 4258 (open circle), and NGC 4596 {8B, filled circle), in the middle panel 
{8B fit) are NGC 1068, NGC 2778, NGC 4258, NGC 4486 {15E, Cap8, filled star) and IC 
U59{15E, filled star), and in the bottom panel {23gal fit) are NGC 1068, NGC 2778, and 
NGC 4258. 
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Ellipticals: Mbh-*(Ry8) Spiral/Lent Bulges: Mbh-*(Rya) All Galaxies: Mbh-*(Rj,/a) 
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Fig. 10. — Same as Figure [5l but for Gravitational Potential at Re/8. 
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Fig. 11. — Same as Figure [9l but for Gravitational Potential at Re/8. The galaxies with 
\residuals\ > 0.6 in all panels are NGC 1068 (open circle), NGC 3384 (open circle) and NGC 
4742 (open circle). Additionally, in the top panel {15E fit): NGC 1023 {8B, filled circle), 
NGC 2787 {8B, filled circle), NGC 4342 (open circle) and NGC 4596 {8B, filled circle), in 
the middle panel {SB fit): NGC 221 {15E,Cap8, filled star), NGC 4486 {15E,Cap8, filled 
star), NGC 4649 {15E, filled star), and IC 1459 {15E, filled star) and in the bottom panel 
{23gal fit): NGC 1023 {8B, filled circle) and NGC 4596 {8B, filled circle) are outliers. It 
is notable that for the bulge {8B) fit the high-mass and low- mass elliptical galaxies are all 
outliers, while for the elliptical fit {15E) fit 3 out of the 8 fitted bulges are outliers; this 
further illustrates the discrepancy between these two sample fits. 
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Fig. 12.— Same as Figure El but for M(10i?e)- 
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Fig. 13. — Same as Figure [9], but for M(lO-Re)- The galaxies with \residuals\ > 0.6 in all 
panels are NGC 821 {15E,Cap8, filled star), NGC 1068 (open circle), NGC 2778 (open star), 
NGC 4342 (open circle) and NGC 4258 (open circle). Additionally, in the top panel {15E 
fit) NGC 7457 {8B, filled circle) is an outlier. 
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Fig. 14. — Results on the calculated gravitational binding energy when the inclination angle 
(top) and observed axis ratio (bottom) are varied over the range of physically allowed values 
for three galaxies: NGC 221 (left), NGC 6251 (middle) and NGC 4486 (right). 



